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ABSTRACT 


The  Physics  Department  at  the  Naval  Postgraduate  School  is  developing  a 
concept  to  overcome  the  problems  that  keep  present  rail  guns  from  being  practical 
weapons.  The  rails  must  be  replaced  often  if  the  rail  gun  operation  is  to  be  continuous. 
Replacing  the  rails  in  present  rail  gun  configurations  is  time  consuming.  The  Physics 
Department’s  design  concept  uses  a  rectangular  barrel  as  part  of  the  solution  to  the 
problem  of  replacing  the  rails.  The  projectile  will  require  flat  surfaces  to  maintain 
electrical  contact  with  the  flat  rails  and  aerodynamic  stabilization  because  of  the  lack  of 
angular  momentum.  This  thesis  develops  one  possible  moctel  of  a  projectile  for  a 
rectangular  barreled  rail  gun,  which  could  be  used  to  replace  the  standard  five-inch  gun 
found  on  most  warships.  The  proposed  projectile  is  successfully  modeled  as  a  five  inch 
projectile  with  flat  areas  planed  onto  opposite  sides  and  long  chord,  short  span  fins 
attached  in  a  cruciform  configuration.  The  computer  programs  used  to  develop  the 
projectile  model  are  included  to  allow  evaluation  of  alternate  configurations. 
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1.  INTRODUCTION 


This  thesis  models  an  aerodynamically  stable  projectile  for  a  rectangular  barreled 
rail  gun.  The  three  con5)uter  programs  used  in  modeling  this  projectile  are  provided  so 
that  the  models  can  be  updated  to  meet  changing  requirements. 

The  first  section  of  Chapters  n  and  IE  present  the  reasonmg  and  assumptions  for 
the  selections  made  from  the  available  alternatives  regarding  that  chapter's  subject.  The 
second  section  of  Ch£q)ters  n  and  III  develop  the  equations  used  in  the  conqjuter  model 
associated  with  that  chapter.  Chapter  IV  uses  the  assumptions  from  the  previous  two 
chapters  gnH  begins  with  four  sections  of  equation  development.  The  final  section  of 
each  chapter  analyzes  the  output  of  the  computer  program  developed  in  the  previous 
sections  from  that  chapter. 

This  thesis  is  multidisciplinary,  with  advisors  in  both  the  Aeronautical 
Engineering  Department  and  in  the  Physics  Department.  Readers  from  different 
disciplines  are  reminded  that  they  may  find  more  detail  relating  to  their  particular  subject 
area  than  they  feel  necessary,  but  readers  from  other  disciplines  may  need  that  detail  to 
understand  what  is  being  done. 

A.  BACKGROUND 

A  rail  gun  uses  the  interaction  between  electrical  currents  and  magnetic  fields  to 
accelerate  a  projectile.  The  rails  in  a  rail  gun  are  electrically  charged  with  the  projectile 
completing  the  electrical  circuit  between  them.  An  electric  current  flows  through  the 
rails  and  projectile  creating  a  magnetic  field.  The  magnetic  field  interacts  with  the 
flowing  electrons,  which  accelerates  the  projectile  out  of  the  barrel.  The  acceleration  is 
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increased  as  the  magnetic  field  is  increased.  The  magnetic  field  increase  as  the  electric 
current  increases. 

The  Physics  Department  at  the  Naval  Postgraduate  School  has  a  rail  quick  change 
concept  for  a  rail  gim.  The  large  electrical  current  used  in  rail  guns  to  achieve  the  large 
accelerations  to  attain  large  muzzle  velocities  in  present  rail  gun  designs  also  vaporizes 
some  of  the  surface  of  the  projectile.  The  vaporized  material  settles  on  the  rails,  fouling 
the  barrel  in  about  five  firings  or  less.  The  rails  must  be  removed  for  cleaning  before  the 
rail  gun  can  be  fired  again.  Present  rail  gun  designs  take  thirty  minutes  of  longer  to 
remove  and  replace  the  rails,  assuming  clean  replacement  rails  are  available. 
Conventionally  fired  guns  can  maintain  a  rate  of  fire  of  several  shots  per  minute  for 
extended  periods.  A  thirty-minute  cease-fire  between  five  shot  salvoes  is  not  an 
acceptable  rate  of  fire  for  a  rail  gun  weapon  system. 

The  Naval  Postgraduate  School  Physics  Department’s  rail  gun  design  concept, 
among  other  things,  will  shorten  the  time  required  to  replace  the  rails  and  thus  increase 
the  firing  rate  in  comparison  with  existing  rail  gun  designs.  The  computer  models 
developed  in  this  thesis  provide  graphs  for  the  analysis  of  the  aerodynamic  effects 
imposed  upon  a  projectile  by  this  alternate  design  as  physical  dimensions  are  changed. 
The  resultant  projectile  performance  can  also  be  analyzed.  Only  those  details  of  the  rail 
gun  design  effecting  the  design  of  the  projectile,  such  as  barrel  shape,  will  be  discussed 
here.  Rail  gun  design  is  such  a  broad  field  of  study  that  all  other  details  of  the  concept, 
such  as  power  requirements  and  building  techniques.  Are  left  to  other  investigators. 

The  most  important  feature  of  the  new  rail  gun  concept  is  its  rectangular  barrel. 
The  projectile,  in  common  with  all  rail  gun  projectiles,  requires  an  electrical  contact 
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suirfece  capable  of  withstanding  the  immense  electrical  currents  used  in  rail  guns.  The 
magnitude  of  these  electrical  currents  depends  on  the  weight  of  the  projectile,  the  desired 
muzzle  velocity,  and  the  length  of  the  barrel.  The  electrical  contact  surfeces  of  the  new 
design's  rectangular  barrel  will  be  flat  rails  on  oppostte  sides  of  the  barrel  Therefore,  to 
maintain  good  electrical  contact  the  electrical  contact  surfece  on  the  projectile  will  have 
to  be  flat.  The  rectangular  barrel,  having  no  rifling,  will  not  impart  the  rotation  that 
provides  the  angular  momentum  that  usually  keeps  projectiles  stable  in  flight.  The  lack 
of  angular  momentum  adds  the  requirement  of  some  form  of  aerodynamic  stabilization 
for  the  projectile. 

B.  BASIC  THEORY 

A  variety  of  aerodynamic  theories,  presented  by  S.  S.  Chin  [Ref.  1]  and  by 
William  Pitts,  Jack  Nielsen,  and  George  Kaattari  [Ref.  2],  were  used  in  development  of 
tihe  conqjuter  models  used  herein  to  compare  projectile  shapes  and  related  performance. 
Slender  body  and  nose/body  theories  from  Reference  (1)  were  used  to  determine  the  skin 
friction  drag,  form  drag,  boattail  drag,  and  base  drag  of  the  body  of  the  projectile.  Airfoil 
theory,  also  from  Reference  (1),  was  used  for  analyzing  small  fins  added  to  the  projectile. 
Figure  (1)  on  the  next  page  shows  these  stabilizing  fins  attached  to  the  projectile  body. 
These  fins  conq)ensate  for  the  lack  of  spin  stabilization  normally  provided  by  the  angular 
momentum  caused  by  the  rotation  of  the  projectile.  The  fins,  however,  also  add  drag. 
Wing/body  interference  fectors  from  Reference  (2)  were  used  in  the  equations  that 
combine  the  aerodynamic  characteristics  of  the  projectile  body  and  fins  to  provide  a  total 
drag  coefficient  for  the  equations  of  motion  and  the  aerodynamic  coefficients  for  the 
stability  study  of  the  complete  projectile. 
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Figure  (1)  Projectile  Schematic 

3.6  2.4 

|<-m  -4 - 18.1  in - — ^.34 


Numerical  integration  theory  is  required  for  the  computer  model  that  calculates 
projectile  trajectory  and  impact  energy.  As  a  projectile  travels  through  space,  many 
factors  interact.  At  any  point  in  time  the  velocity  is  changing  due  to  the  effects  of  drag, 
which  in  turn  changes  the  distance  the  projectile  will  have  traveled  at  the  next  point  in 
time  and  the  velocity  it  has  when  it  gets  there.  To  further  complicate  matters,  the  drag 
itself  changes  with  the  velocity  and  with  the  air  density,  which  is  changing  with  the 
vertical  position  of  the  projectile.  With  so  many  interactions  between  variables  a  closed 
solution  to  the  integration  of  the  equations  of  motion  is  difficult,  consequently  numerical 
integration  schemes  are  required  to  solve  the  equations  of  motion.  The  Euler-Cromer 
method  was  selected  for  the  numerical  integration  of  the  equations  of  motion  because  of 
its  stability  in  long  calculations.  The  Euler-Cromer  method,  occasionally  referred  to  as 
the  last  point  method,  uses  a  loop,  shown  in  simplified  form  in  Figure  (2)  as  a  flow  chart. 
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The  program  loop  takes  the  final  position  of  the  projectile  firom  the  previous  time  step 
and  calculates  the  velocity  at  that  point  and  then  uses  that  velocity  to  calculate  the  next 
point  in  the  projectile's  trajectory.  The  accuracy  of  the  integration  for  most  numerical 
integration  methods  is  increased  as  the  time  steps  are  made  smaller,  but  the  round-ofif 
error  often  increases  with  each  calculation.  The  Euler-Cromer  method  is  stable  in  long 
calculations  because  the  round-off  errors  tend  to  cancel  each  other  so  the  calculation  does 
not  diverge. 

Figure  (2)  Euler  Cromer  Numerical  Xntegration  Flow  Chart 
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n.  SIZING 


There  are  many  different  ways  of  fitting  a  projectile  in  a  rectangular  barrel.  At 
the  time  of  this  writing  the  design  of  the  rail  gun  is  still  very  conceptual,  thus  there  are 
few  physical  limits  placed  upon  the  projectile. 

A  MATLAB  conqjuter  prograna,  provided  in  Appendix  (A),  is  used  to  model  the 
different  possible  projectile  nose/body  configurations  as  the  radius  of  the  projectile  is 
changed.  This  MATLAB  computer  program  will  be  referred  to  as  the  configuration 
conparison  program.  Since  some  of  the  design  alternatives  chosen  may  become  invalid 
as  the  rail  gun  design  is  formalized,  the  equations  for  the  configuration  conparison 
program  have  been  generalized.  Future  users  of  the  program  will  be  able  to  change 
inputs,  such  as  the  nose  fineness  ratio  or  the  explosive  payload  of  the  projectile  model,  to 
suit  new  requirements  and  the  model  will  still  be  correct. 

A.  CONFIGURATION  ALTERNATIVES 

The  configuration  choices  made  from  among  the  different  projectile  design 
options,  e.g.,  ogive  verses  spherical  nose,  body  diameter,  etc. . .  will  determine  the 
projectile's  physical  dimensions.  The  physical  dimensions,  such  as  cross  sectional  area, 
will  have  a  direct  effect  on  the  aerodynamic  characteristics  of  the  projectile. 

1.  Electrical  Contact  Surface 

The  projectile  requires  flat  electrical  contact  surfeces  while  in  the  rail  gun  barrel. 
Two  methods  to  do  this  are  immediately  evident;  place  the  projectile  in  a  flat-sided  sabot 
or  put  flat  sides  on  the  projectile.  Flat  sides  on  the  projectile  for  the  electrical  contact 
area  are  used  in  the  model  because  of  the  two  major  problems  with  sabots  outlined  below. 
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A  sabot  is  a  spacer  placed  around  a  projectile  to  allow  the  projectile  to  be  fired 
from  a  larger  or  differently  shaped  barrel.  Sabots  are  usually  light  and  often  used  in 
conventional  guns  when  using  a  larger  powder  charge  than  would  otherwise  be  possible 
behind  a  smaller  projectile  to  increase  the  projectile's  muzzle  velocity.  A  sabot  seems 
attractive  from  an  aerodynamics  consideration  because  the  projectile  can  remain  truly 
cylindrical.  A  cylindrical  projectile  should  have  reduced  airflow  separation  and  surface 
flow  field  turbulence.  The  sabot  in  the  rail  gun,  however,  has  to  conduct  the  electrical 
current  required  to  accelerate  the  projectile  without  melting.  This  will  require  a 
substantial  piece  of  metal,  which  equates  to  a  large  mass  and  thus  a  large  inertia.  This 
large  mass  will  have  to  be  accelerated  with  the  projectile  and  the  large  inertia  means  it 
will  take  time  to  stop  once  it  has  left  the  barrel 

Accelerating  the  sabot  becomes  an  energy  consideration.  The  energy  at  the  target 
is  half  the  impact  mass  times  the  impact  velocity  squared.  Since  a  sabot  separates  from 
the  projectile  immediately  after  leaving  the  barrel  its  mass  is  not  present  at  projectile 
impact  so  the  energy  used  to  accelerate  the  sabot  is  wasted. 

The  large  inertia  associated  with  the  sabot  is  a  problem  because  the  projectile  will 
separate  from  the  sabot  as  they  both  exit  the  barrel  of  the  rail  gun  and  the  sabot  will  land 
somewhere  down  range.  The  light  plastic  sabots  in  the  Phalanx  gun,  for  example,  have 
comparatively  low  velocity,  high  drag,  and  low  inertia  so  they  travel  only  a  short  distance 
and  fall  harmlessly.  The  sabot  from  the  proposed  rail  gun's  projectile  will  have  high 
drag,  but  it  will  have  high  velocity  and  high  inertia.  This  means  that  on  long  range  shots 
the  sabot  may  travel  several  miles  down  range  and  land  with  a  high  impact  velocity. 
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Depending  on  the  location  of  Siendly  forces  and  neutral  populations  along  the  line  of 
&e,  where  the  sabot  comes  down  could  be  as  important  as  where  the  projectile  impacts. 
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Figure  (3)  Nose  and  Body  Configurations 
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Two  methods  for  providing  flat  electrical  contact  areas  on  the  projectile  are 
considered.  F^ure  (3)  above  illustrates  both  of  these  methods.  The  first  method  is  to 
slice  a  cylindrical  projectile  lengthwise  down  the  middle  and  insert  a  rectangular  wedge 
of  ^ropriate  width  with  the  fix)nt  of  the  wedge  matching  the  curvature  of  the  nose.  The 
second  option  for  providing  the  electrical  contact  surfece  is  to  shave  down  opposite  sides 
of  a  cylindrical  projectile  until  a  sufficiently  wide  fiat  area  is  planed  out.  This  reduces 
cross  sectional  area,  but  lengthens  the  projectile  for  the  same  ejqjlosive  payload.  Neither 
has  an  intuitive  advantage  over  the  other  so  both  are  modeled  in  the  configuration 
con^arison  program  and  a  choice  will  be  made  by  a  trade  study  of  the  output. 
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The  width  of  the  electrical  contact  area  must  be  considered.  Making  the  electrical 
contact  surfece  smaller  makes  the  body  more  cylindrical,  which  should  reduce  the 
separation  and  turbulence  of  the  airflow  around  the  body  and  therefore  has  aerodynamic 
benefit.  A  smaller  contact  area,  however,  means  that  the  electric  current  has  to  go 
through  a  smaller  area,  which  at  the  high  currents  required  for  rail  gun  operation  leads  to 
vaporization  and  spalling  of  the  projectile’s  flat  surfece.  Making  the  contact  surfece 
wider  has  the  opposite  effects  of  increasii^  separation  and  turbulence  of  airflow  aroimd 
the  body  and  reducing  current  density.  The  two-inch  wide  contact  area  shown  in  Figure 
(3)  is  selected  somewhat  arbitrarily  because  it  provides  a  sufficiently  large  electric 
contact  area  and  the  projectile  remains  fairly  cylindrical.  The  width  of  the  contact  area  is 
one  of  the  physical  parameters  that  can  be  changed  by  the  user  in  the  input  section  of  the 
configmation  conqjarison  program  to  account  for  changing  requirements.  Since  the 
contact  surfece  extends  the  entire  length  of  the  body  a  narrow  contact  area  may  still 
provide  a  large  electrical  contact  surfece  on  both  sides  of  the  projectile.  The  drag 
equations  in  the  configuration  conq>arison  program  are  written  always  assuming  turbulent 
flow.  Turbulent  flow  provides  the  worst  drag  which  somewhat  masks  the  effect  of  the 
flattened  region.  As  specific  parameters  of  the  rail  gun  design  are  set  and  specific 
materials  are  chosen  a  smaller  contact  area  may  be  possible. 

2.  Nose  Type 

There  are  several  possible  nose  types  for  the  projectile.  The  proposed  maximum 
muzzle  velocity  of  the  rail  gun  design  is  Mach  8.9.  The  nose  will  experience  extreme 
temperatures  for  the  first  several  seconds  after  launch  because  of  skin  fiiction  heating  in 
the  lower  atmosphere  and  the  shock  front  heating  the  air  around  the  projectile.  The  high 
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temperatures  will  last  until  the  projectile  velocity  is  reduced,  thereby  reducing  the 
strength  of  the  shock,  or  the  projectile  reaches  the  upper  atmosphere,  where  the  lower  air 
density  will  reduce  friction. 

N.  Kemp  and  F.  Riddell  [Ref.  3]  show  the  heating  of  a  projectile’s  nose  varies 
inversely  with  the  square  root  of  the  radius  of  the  nose’s  point  as  seen  in  equation  (2.1). 

Q  =  20800(R„)-‘^'(p/po)' (2.1) 
The  other  factors  in  the  equation  are  environmental  factors  independent  of  the  projectile. 
The  temperature  change  of  the  nose  is  determined  by  equation  (2.2). 

Mcp(dT/dt)  =  AQ  (2.2) 

The  factor  AQ  is  the  amoimt  of  heat  being  added  by  skin  friction  and  the  shock  front 
minus  the  heat  being  radiated  away  by  the  nose  to  the  atmosphere  or  conducted  away 
from  the  nose  to  other  parts  of  the  projectile.  A  very  pointed  nose,  meaning  a  small  nose 
point  radius,  provides  less  drag,  will  experience  more  heating,  and,  because  it  has  less 
mass  over  which  to  distribute  that  heat,  risks  melting  and  deforming.  A  very  blunt  nose, 
meaning  a  large  nose  point  radius,  has  higher  form  drag,  will  experience  less  heat,  and, 
because  it  has  more  mass  over  which  to  distribute  that  heat,  is  at  reduced  risk  of  melting 
and  deformation.  The  rate  at  which  the  nose  is  heated  and  Ihe  thermal  conductivity  and 
emissivity  of  its  constituent  material  determine  the  amount  of  mass  required  to  absorb  the 
heat  without  melting.  If  the  nose  is  heated  faster  than  it  can  radiate  or  diffuse  the  heat 
away  if  may  melt  or  even  vaporize. 

A  tangent  ogive  with  a  fineness  ratio,  nose  length  divided  by  the  maximum 

f 

diameter,  of  unity  is  used  in  the  configuration  comparison  computer  program.  The 
fineness  ratio  is  another  of  the  physical  parameters  that  the  user  can  change  in  the  input 
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section  of  the  configuration  con^arison  program.  The  user  can  select  a  fineness  ratio  to 
minimize  the  drag  while  maintaining  heat  capabilities  for  the  conditions  ejq)ected.  A 
spherical  nose,  alternative  to  the  tangent  ogive,  can  also  be  modeled  for  comparison. 
Figure  (3)  on  a  previous  page  illustrates  the  two  methods  of  providing  the  two-inch  wide 
contact  sur&ce  with  both  nose  types  and  a  two  and  one  half-inch  radius  of  body 
curvature. 

3.  Component  Material 

The  material  of  which  the  projectile  is  made  has  important  consequences  on  the 
eventual  performance  of  the  projectile.  The  material  has  to  be  a  good  conductor  of  both 
electricity  and  heat.  A  high-density  material  offers  several  advantages  resulting  in  longer 
ranges  and  higher  impact  velocities. 

The  biggest  benefit  of  using  high  density  material  is  the  increased  mass  per  unit 
volume.  The  drag  force  is  determined  by  the  geometry  of  the  projectile  and  the 
environment  in  which  the  projectile  is  placed.  The  force  on  the  projectile  equals  its  mass 
times  its  acceleration.  If  the  geometry  of  the  projectile  is  not  changed  so  the  drag  force  is 
constant  at  a  point  in  space,  but  the  mass  is  increased  then  the  acceleration,  deceleration 
in  the  case  of  a  drag  force,  must  decrease  as  shown  in  Figure  (4)  on  the  next  page. 
Therefore  the  effect  of  the  drag  force  on  the  velocity  of  the  projectile  is  inversely 
proportional  to  the  mass  of  the  projectile.  The  terminal  velocity  of  the  projectile  is 
directly  proportional  to  the  mass  of  the  projectile  as  calculated  by  equation  (2.3). 

Vt  =  (2mg/pScsCi>)*^  (2.3) 
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p  is  the  ambient  air  density,Scs  is  the  projectiles  cross  sectional  area,  and  Cd  is  the  total 
drag  coefficient.  Another  benefit  of  high-density  material  is  the  increased  projectile 
penetration  capabilities  against  hard  targets. 
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Tungsten  is  very  dense  and  is  already  used  for  armor  piercing  projectiles. 
Tungsten  is  a  good  choice  for  this  projectile  and  is  used  in  the  configuration  con^jarison 
program  for  these  reasons.  The  density  of  the  component  material  can  be  selected  by  the 
user  in  the  input  section  of  the  program. 

As  with  most  thing  there  are  always  disadvantages  to  any  option.  For  the 
reason  the  increased  mass  of  a  projectile  made  of  tungsten  resists  the  deceleration  caused 
by  drag,  it  also  resists  the  acceleration  by  the  electric  and  magnetic  fields  in  the  rail  gun, 
requiring  higher  currents  and  longer  barrels.  Shipboard  portability  is  another  major  issue 


for  very  heavy  projectiles  as  the  individual  projectiles  may  no  longer  be  handled  by  a 
single  sailor. 

The  disadvantages  of  large  mass  will  probably  drive  an  actual  projectile  to  be 
made  of  a  lighter  alloy,  possibly  keeping  an  all  tungsten  nose  for  its  penetration 
capabilities.  There  are  limits,  however,  on  how  light  a  projectile  can  be  made  for  a 
constant  volume.  As  seen  in  Figure  (4)  the  deceleration  caused  by  drag  will  be  greater 
for  a  lighter  projectile  and  therefore  the  range  and  impact  velocity  will  be  little  better  than 
a  conventionally  fired  projectile. 

The  configuration  comparison  program  can  use  the  density  of  any  material 
without  efifecting  the  accuracy  of  the  modeling.  The  equations  in  the  con:q)uter  program 
that  models  the  trajectory  of  the  projectile,  provided  in  Appendix  (B)  and  fully  developed 
in  chapter  four,  are  written  assuming  that  the  projectile  remains  supersonic  all  the  way  to 
impact.  If  the  projectile  becomes  subsonic  anywhere  in  the  trajectory  the  calculations  in 
the  trajectory  program  diverge  because  several  of  the  formulas  are  not  defined  for  the 
subsonic  regime.  Therefore  the  projectile's  component  material  must  be  sufficiently 
dense  to  keep  the  projectile  supersonic  until  impact. 

4.  Explosives 

A  solid  metal  projectile  at  high  velocities  has  a  great  deal  of  kinetic  energy,  but  at 
long  ranges  where  the  projectile  has  slowed  it  may  need  to  carry  ejq)losives  to  be 
effective.  Against  hardened  targets,  such  as  ships  or  concrete  structures,  the  sudden 
deceleration  of  the  projectile  at  impact  transfers  that  kinetic  energy  into  the  target. 

Against  soft  targets,  such  as  wooden  buildings  or  troops,  there  is  no  mechanism  to  rapidly 
convert  the  energy  of  the  projectile  into  destructive  power  at  the  target.  The  destructive 
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power  is  the  impact  energy  divided  by  the  time  over  which  it  is  transferred.  The  shock 
wave  of  the  projectile,  if  it  is  still  at  a  high  Mach  number  when  it  reaches  the  target,  may 
some  damage  because  of  the  pressure  increase  with  the  passing  of  the  shock  wave. 
The  higher  the  Mach  number  of  the  projectile  the  stronger  the  shock  will  be.  Figure  (5) 
on  the  next  page,  computed  with  an  equation  presented  by  R.  Zucker  [Ref.  4]  used  to 
formulate  the  normal  shock  tables,  shows  the  upper  bound,  since  normal  shocks  are  the 
strongest  sbfx^ks^  of  the  pressure  ratio  as  the  Mach  number  is  increased.  The  pressure  rise 
caused  by  the  shock  wave  quickly  dissipates  with  distance  from  the  projectile  and  the 
obliquity  of  the  shock  front.  The  projectile  itself  will  only  punch  a  small  hole  through  the 
walls  of  a  wooden  building  or  into  the  soft  dirt  like  a  present  day  armor  piercing  shell. 

Against  eitha:  hard  or  soft  targets  at  long  ranges  the  projectile  will  decelerate  to 
terminal  velocity  before  impact  and  may  not  have  much  more  kinetic  energy  than  a 
convetional  projectile.  If  the  rail  gun  is  to  be  a  practical  we^on  it  must  be  able  to 
perform  at  least  as  well  as  the  current  weapons,  such  as  the  five  inch  gun  it  may  someday 
replace.  This  means  the  proposed  projectile  will  have  to  carry  an  e3q)losive  payload 
similar  to  a  convetional  projectile.  Twelve  pounds  have  been  selected  as  the  payload 
mass  for  the  model  because  that  is  the  mass  of  explosive  material  in  a  five  inch  shell 
The  e3q>losive  nnmpmmd  C4  is  used  in  the  conq)uter  model  because  the  density  of  C4  is 
about  average  for  modem  high  energy  ejqplosives.  The  explosive  payload  and  the  density 
of  the  selected  explosive  are  both  variables  which  may  be  changed  by  the  user  in  the 
input  section  of  the  computer  model  The  length  of  the  projectile  is  calculated  from  the 
volume  required  to  contain  the  selected  mass.  Substituting  a  different  e:q)losive 
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Figure  (5)  PrassureRaflolbr  Passage  of  Normal  Shock 


compound  of  similar  density  at  the  same  explosive  payload  should  only  make  a  small 
difference  in  the  length  of  the  projectile. 

5.  Casing  Thickness 

The  explosive  for  the  projectile  must  obviously  be  encased  within  the  projectile  so 
the  question  of  casing  thickness  arises.  A  casing  thicker  than  a  convetional  projectile 
would  provide  more  material  to  help  absorb  the  e^qiected  aerodynamic  heat  transfer.  The 
projectile's  cross  sectional  area  and  surface  area,  however,  will  be  increased  for  a 
constant  explosive  payload  which  increases  the  drag.  The  case  thickness  is  yet  another 
variable  that  a  user  can  select  in  the  model's  input  section.  A  half  inch  casing  thickness 
was  used  in  the  model  presented  in  the  final  section  of  this  chapter. 


Additional  options  regarding  materials  emerge  at  this  point.  Using  a  single 
material  for  the  entire  projectile  eases  manufacture  since  the  nose  and  body  can  be  cast  as 
a  single  unit,  so  the  model  uses  a  single  material,  but  the  conputer  model  can  easily  be 
altered  to  make  the  nose ,  end  plate  and  the  casing  out  of  different  materials.  This  is 
significant  because  the  nose  is  not  in  contact  with  the  conducting  rails  rails  so  its  electric 
conductivity  isn't  very  important  and  the  casing  will  be  the  main  conductor  of  electricity 
in  the  barrel  so  its  conductivity  is  of  great  importance. 

B.  DIMENSIONS  AND  FORGES 

This  section  details  the  development  of  the  equations  for  the  configuration 
comparison  program  that  define  the  physical  dimensions  of  the  projectile's  nose  and 
body.  Many  of  the  calculations  are  complicated  by  keeping  the  physical  dimensions  as 
variables  to  allow  the  flexibility  in  the  computer  model  discussed  in  the  previous  section. 

1.  Nose 

The  mass  and  surface  area  of  the  nose  are  the  important  parameters  for 
determining  the  the  effects  of  the  aerodynamic  properties  in  the  equations  of  motion. 
Refering  back  to  Figure  (3)  will  help  the  reader  visualize  the  nose  configurations.  To 
determine  the  mass,  the  volume  of  the  nose  and  the  mass  density  of  its  conqjonent 
material  are  required.  Once  the  volume  is  determined  the  mass  can  be  calculated  by 
multiplying  the  volume  and  the  density. 

The  geometry  of  the  wedge  insert  configuration  is  simpler,  making  the  required 
parameters  easier  to  calculate.  The  volume  of  a  half  sphere  and  conical  ogive  are  known 
formulas,  the  formula  for  the  ogive  can  be  found  in  Appendix  A  of  Reference  (1).  The 
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inserted  wedge  for  the  spherical  nose  is  just  a  half  cylinder  of  the  same  radius  and  desired 
width.  Using  r  as  the  radius  of  the  spherical  nose  and  h  as  the  width  of  the  contact  area  in 
equation  (2.4). 

NVs^v  =  (2/3)7rr^  +  (l/2)7tI^h  (2.4) 

The  calculation  for  the  ogive  wedge  configuration  is  a  little  more  involved.  Calculus 
provides  a  convenient  solution  for  the  area  of  a  section  of  a  circle,  shown  in  Figure  (6) 
below.  To  get  the  volume  we  multiply  the  planform  area  that  matches  the  nose  by  two 
and  the  desired  width.  Equation  (2.5)  uses  L  as  the  nose  length,  R  as  the  radius  of 
curvature  of  the  nose,  and  r  and  h  as  above. 

NVo/w  =  7r[L(2R2.2Rr+r2)  -  lVS  -  (R-r)L(R2-L2)‘^  -  (R.-T)Rhm^  (L/R)]  + 

\R\n/2-Q)  -  (R-r)^cote]h  -  (2.5) 

Figure  (6)  Calculating  Radius  of  Curvature  and  Planform  Area  for  Ogive 

r  =  Radtus  of  ProjectQe 
L  =  Nose  Length 


solving  the  integrals 


|lR^(^-e)-(R-r)2cot0] 

The  volume  of  the  shaved  side  projectile  configuration  is  easily  understood  in 
concept,  but  mathematicaly  complicated  to  confute.  The  concept  is  simply  to  determine 
the  volume  for  the  complete  nose  shape  and  subtract  the  volume  of  the  piece  that  is 
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shaved  away.  The  removed  volume  is  determined  by  solving  the  tr^le  integral  in 
equation  (2.6). 


h/2  7t/2  R 

RV=  1  J  J  r'dr'de’dz  (2.6) 

0  cos'‘(fe/R)  (R^-i?z^)^'^csc0' 

The  term  f  is  the  nose  fineness  ratio.  Completing  the  integration  with  variables  for  the 

limits  of  integration  requires  an  approximation.  Simpson's  Rule  is  used  to  get  the 

solution  of  the  integral  above  for  equation  (2.7): 

RV  =  rW  -  (R2h/4)cos  ^(fli/2R)  -  (RV2f)(l-(fh/2R)^)^^  +  RV2f  -  (2.7) 

(Rfh/48)[(h/2)csc(cos'Vh£'8R))  +  (h/2)csc(cos-^(hf4R))  + 

(3h/2)csc(cos'\3h£'8R))  +  (h/2)csc(cos‘Vhf'2R))]  + 
(^h/48R)[(h^/128)csc(cos-‘(hfi'8R))  +  (h^/32)csc(cos'*(hf'4R))  + 

(27hVl28)csc(cos-^(3hj6'8R))  +  (hV8)csc(cos‘(h£'2R)) 


Multiplying  this  volume  by  four  gives  the  total  volume  lost  from  the  nose  by  planing  the 
flat  sides  on  to  the  projectile.  A  spherical  nose  is  the  limiting  case  of  the  ogive  nose 
where  the  fineness  ratio  equals  one  half,  the  radius  of  curvature  of  the  nose  and  the  length 
of  the  nose  equal  the  radius  of  the  projectile.  The  volume  equation  for  the  shaved  side 


spherical  nose  configuration  is  equation  (2.8). 

NVs/s  =  (2/3)to^  -  R"lm/2  +  R''hcos-‘(h/4R)  +  4R^(l-(h/4R)")'"'  -  R^  +  (2.8) 

(Rh/24)[(h/2)csc(cos'^(h/16R))  +  (h/2)csc(cos'^(h/8R))  + 

(3h/2)csc(cos'*(3h/16R))  +  (h/2)csc(cos'^(h/4R))]  - 
(h/96R)[(hVl28)csc(cos  *(h/16R))  +  (h^/32)csc(cos-*(h/8R))  + 

(27h^/128)csc(cos'T3h/16R))  +  (h^/8)csc(cos’‘(h/4R)) 

The  volume  for  the  shaved  side  ogive  nose  configuration  is  shown  in  eqxiation  (2.9); 

NVo/s  =  7t[L(2R^-2Rr+T2)  -  V/3  -  (R-r)L(R^-L2)^^  -  (R-r)sm‘(L/R)]  -  (2.9) 

R^lm/2  +  R2hcos'^(fh/2R)  +  (2R^/f)(l-(fh/2R)^)^^  -  2RVf + 
(Rfh/12)[(h/2)csc(cos'*(h£'8R))  +  (h/2)csc(cos"*(hfi'4R))  + 
(3h/2)csc(cos-^(3hfr8R))  +  (h/2)csc(cos ‘(hfr2R))]  - 
(f®h/12R)[(h^/128)csc(cos-*(hfr8R))  +  (h^/32)csc(cos-‘(hfr4R))  + 
(27h^/128)csc(cos-H3hfi'8R))  +  (h2/8)csc(cos'*(hfr2R)) 
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The  surface  area  of  the  various  nose  configurations,  needed  to  determine  skin 
friction  drag.  Once  again  the  wedge  insert  configuration  surface  areas  were  easier  to 
calculate.  The  spherical  nose  with  wedge  insert  is  the  surface  area  of  a  half  sphere  plus 
the  surface  area  of  half  a  cylinder  of  matching  radius  and  desired  electrical  contact  area 
width.  The  formula  for  the  spherical  nose  with  wedge  insert  is  given  in  equation  (2.10): 

Sns/w  =  27tr^  +  Tirh  (2.10) 

The  equation  for  the  surfece  area  of  the  ogive  nose  is  a  surface  of  revolution  and  is 
determined  by  solving  the  integral  in  equation  (2.1 1): 

2fr 

So  =  1  27t[(RV)*^  -  (R-r)][l  V/(RV)]’'^  dy  =  27t[2frR  -  (2.11) 

0  R(R-r)sm*(2fr/R)] 

The  surface  area  of  the  inserted  wedge  is  the  desired  electrical  contact  area  width  and  the 

arc  length,  \^ch  is  determined  from  the  radius  of  curvature  of  the  nose  and  6  as 

calculated  in  Figure  (6),  doubled  to  accoimt  for  the  upper  and  lower  surface  of  the  wedge. 

The  surface  area  of  the  ogive  nose  with  wedge  insert  is  the  sur&ce  area  of  the  ogive  nose 

plus  the  surface  area  of  the  wedge  insert,  shown  in  equation  (2.12). 

S„o/w  =  27r[2frR  -  R(R-r)sm  *(2fr/R)]  +  2hRe  -  (2. 12) 

The  surface  area  of  the  shaved  side  configuration  is  calculated  in  much  the  same 
manner  as  the  volrune  was  determined.  The  sur&ce  area  of  the  removed  section  is 
subtracted  form  the  surface  area  of  the  whole  surfece,  but  now  the  surface  area  of  the 
remaining  flat  surface  must  be  added  back  to  the  surface  area.  For  the  shaved  side 
spherical  nose  one  can  use  the  symmetry  as  we  did  in  the  volume  calculation.  The 
removed  surface  area  can  be  considered  a  surfece  of  revolution  and  determined  by 
solving  the  integral  in  equation  (2. 1 3). 
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.5h 

fo  27ty(l-Kdy/dx)2)^'^  dx  =  27t[{-.5r  lii(.5h  +  (rf.25h2)*^)}  + 

{.25h(r+.25h2y'^}  +  {.5rln(r)}]  (2.13) 

The  integrand  in  equation  (2.13)  is  y  =  (r^-x^)^^  and  the  limits  of  integration  are  from  the 
centerline  of  the  projectile  to  the  edge  of  the  flat  electrical  contact  surfece.  Each  of  the 
resulting  flat  areas  left  on  the  nose  is  a  half  circle  with  a  radius  of  one  half  the  electrical 
contact  area  width.  Collecting  all  the  terms  for  the  surface  of  the  half  sphere,  the  surfttce 
of  the  removed  portion,  and  the  surftice  of  the  remainii^  flat  spots  results  in  equation 
(2.14). 

Sws  =  27tr^  -  27i[{-.5r  ln(.5h+(r+.25h2)‘^)}  +  (2. 14) 

.25h(r+.25h^)‘'^}  +  {.5r  In(r)}]+  7r(h/2)^ 

The  shaved  side  ogive  nose  is  harder  to  calculate  because  of  the  lack  of  symmetry.  The 
removed  surftice  area  for  the  shaved  side  ogive  nose  is  four  times  the  value  of  the  double 
integral  solved  in  equation  (2.15). 

h/2  (7t/2  -  cos'*(fe/r)) 

J  J  R dO dz  =  R[7ih/4 -  (h/2)cos‘(fh/2R) -  (2.15) 

0  0  (R/f)(l-(fh/2R)^)^'^  +  R/fl 

The  flat  area  left  on  the  nose  when  the  sides  ^e  shaved  down  is  once  again  shaped  like 

the  area  in  Figure  (6).  The  surface  area  of  the  ogive  nose  with  shaved  sides  configuration 

is  the  combination  of  the  surfece  area  of  the  ogive  nose  minus  the  surfece  area  of  the 

removed  section  plus  the  surfoce  area  of  the  flat  spots  left  on  the  nose  when  the  shaved 

area  is  remove.  The  ogive  nose  with  shaved  sides  surfece  area  formula  is  given  by 

equation  (2.16). 
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(2.16) 


S„o/s  =  27i[2frR - R(R-r)sm*(2fr/R)]  - R[7ih -  (2h)cos^(fh/2R)  - 
(4R/f)(l-(fh/2R)2)‘^  +  4R/f]  +  2[R^(jt/2-cos-^(h/R))  - 

(R-h/2)^cot(cos'^(h/R))] 

2.  Body 

The  cross  sectional  area,  length,  mass,  and  surface  area  are  the  important  physical 
characteristics  of  the  projectile  body.  The  cross  sectional  area  of  the  wedge  insert 
configuration  is  the  area  of  a  circle  who's  radius  is  the  radius  of  the  projectile  plus  the 
area  of  a  rectangle  who's  height  is  twice  the  radius  and  who's  width  is  the  width  of  the 
contact  area,  as  shown  in  equation  (2.17). 

Scsw  =  7tr^  +  2rh  (2.17) 

The  cross  sectional  area  of  the  shaved  side  configuration  is  the  area  of  a  circle  minus  four 
of  the  planform  areas  that  were  calculated  for  the  shaved  side  spherical  nose,  shown  in 
equation  (2.18). 

Soss  =  -  2[r^(7r/2-0)  -  (r-rO^cotG]  (2.18) 

The  length  of  the  projectile's  body  depends  on  how  much  explosive  is  required. 
First  the  cross  sectional  area  of  the  ejqjlosive  must  be  determined  for  the  desired 
projectile  radius  and  skin  thickness.  The  same  equations  as  were  used  for  the  body's 
cross  sectional  area  can  be  used  to  determine  the  ejqjlosive's  cross  sectional  area  by 
reducing  the  radius  used  in  the  equation  by  the  desired  skin  thickness.  The  length  can 
now  be  determined  by  dividing  the  desired  explosive  payload  mass  by  ejq>losive  cross 
sectional  area  and  the  density  of  the  chosen  explosive.  Equation  (2.19)  uses  me  as  the 
mass  of  the  ejqjlosive,  pe  as  the  density  of  the  explosive,  and  Secs  as  the  appropriate 
e)q)losive  cross  sectional  area  for  the  wedge  or  shaved  configuration. 

Lb  nie/peSecs  (2.19) 
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The  mass  of  the  projectiles  body  is  the  sum  of  the  explosive  mass  plus  the  mass  of 
the  casing  material.  To  determine  the  mass  of  the  casing  it  is  easiest  to  find  the  volume 
of  the  material  by  subtracting  the  explosive  cross  sectional  area  from  the  projectile's  cross 
sectional  area  and  multiplying  by  the  length  of  the  explosive  and  density  of  casing 
Equation  (2.20)  uses  the  variables  from  above  and  Spcs  for  projectile  cross 
sectional  area,  po  for  the  density  of  the  casing  material.  Lb  for  the  body  length,  and  mb  for 
mass  of  the  body. 

nib  ~  nie  +  LbPc(Spcs  ”  Secs)  (2.20) 

Of  course  the  length  and  cross  sectional  areas  must  all  be  for  the  same  wedge  or  shaved 
configuration. 

The  body's  sur&ce  area  for  the  wedge  insert  configuration  is,  again,  easily 
determined.  It  is  the  surfece  area  of  a  cylinder  plus  the  upper  and  lower  contact  surfece 
which  is  the  width  of  the  contact  area  times  the  length  of  the  body  as  shown  in  equation 
(2.21). 

Sbw  =  2OTLbw  +  2hLbw  (2-21) 

The  shaved  side  configuration  requires  reference  back  to  Figure  (3).  In  order  to  get  the 
arc  length  to  determine  the  sur&ce  area  the  figure  should  be  viewed  on  its'  side.  For  this 
case  we  know  that  Rsin0  equals  half  the  contact  area  width  where  R  is  the  radius  of  the 
projectile.  Solving  for  0  and  multiplying  by  the  radius  determines  half  the  arc  length  for 
one  side  of  the  projectile.  Multiplying  the  arc  length  by  four  and  by  length  of  the  shaved 
side  body  gives  the  surfece  area  of  the  curved  portions  of  the  projectile  body.  The  flat 
surfeces  on  the  projectile  are  once  again  the  contact  area  width  times  the  body  length  for 
the  shaved  side  configuration.  Combining  the  surfece  areas  for  the  shaved  side 
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configuration  we  get  equation  (2.22),  the  total  surface  area  for  the  body  in  this 
configuration: 

Sbs  =  4rLbsSm*(h/r)  +  2hLbs  (2.22) 

3.  End  Plate 

The  end  plate  is  assumed  to  be  the  same  thickness  as  the  casing,  although  this  can 
be  changed  in  the  programs  if  desired.  Since  the  thickness  is  set  the  only  parameters  to 
find  are  the  surfece  area  and  mass.  The  surfece  area  can  be  easily  calculated  by  taking 
the  equation  for  the  body  surface  area  for  the  matching  configuration  from  the  previous 
section  and  replacing  the  body  length  with  the  skin  thickness.  The  mass  is  just  as  easily 
determined  by  multiplying  the  cross  sectional  area  of  the  matching  body  configuration  by 
the  skin  thickness  and  component  material  density. 

4.  Drag 

The  form  and.  skin  fiiction  drag  for  the  different  projectile  configurations  are 
computed  using  Reference  (1)  to  compare  the  drag  characteristics  of  each.  The  dynamic 
head  is  required  for  drag  calculations  and  is  computed  in  the  model  using  standard  sea 
level  air  density  and  the  selected  mu2zle  velocity  with  equation  (2.23): 

q  =  (l/2)pv2  (2.23) 

The  skin  friction  drag  is  confuted  using  a  generic  high  velocity  skin  friction  drag 
coefficient  form  Reference  (1)  of  .002.  The  generic  skin  fiiction  drag  coefficient  is  an 
acceptable  approximation  because  the  large  velocity  of  the  projectile  is  the  same  for  all 
cases  and  dominates  the  numerator  of  the  Reynolds  number  calculation,  shown  in 
equation  (2.24)  used  to  determine  an  exact  skin  friction  drag  coefficient,  shown  in 
equation  (2.25). 
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Re  =  pvLt/|i  (2.24) 

Cf  =  .455  [Logio(Re)]-^  ^*  (2.25) 

The  skin  friction  drag,  equation  (2.26),  is  the  skin  friction  drag  coefficient  multiplied  by 


the  total  surface  area  and  the  dynamic  head. 

CDf=CfSq  (2.26) 

The  form  drag  calculation  has  a  few  more  intermediate  steps.  The  form  drag  coefficient 
multiplied  by  the  projectile's  cross  sectional  area,  Spcs,  and  the  djTianiic  head  previously 
calculated  determines  the  form  drag  in  equation  (2.27): 

CDF  =  CF(Spcs)q  (2.27) 

The  form  drag  coefficient  is  calculated  with  equation  (2.28): 

Cf  =  (Ap/q){l-[(196(L„/d)"-16)/(14(M+18)(L„/d)")]}  (2.28) 

The  conical  shock  pressure  coefficient  is  calculated  using  the  enqjirical  formula 
given  in  equation  (2.29). 

Ap/q  =  [.083-K.096/M2)](cy/10)*-^**  •  (2.29) 

The  semivertex  angle,  ct,  of  a  reference  cone  corresponding  to  the  ogive  nose  is 
required  to  determine  the  pressure  rise  across  a  conical  shock  and  is  deteraiined  in 
equation  (2.30)  by  taking  the  arctangent  of  one  divided  by  twice  the  fineness  ratio. 

a  =  tan‘*(l/2Q  (2.30) 

B.  RESULTS 

The  configuration  conqjarison  conputer  program  computes  the  physical  data  for 
the  projectile  over  a  range  of  projectile  radii  selected  by  the  user  for  the  different 
projectile  configurations  and  plots  the  results  for  comparison.  Figures  (7)  through  (12). 
on  the  next  several  pages  presents  these  plots  for  one  con^uter  program  run  with  the 
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Form  Orag/io  *  (kg-m/»^) 


Figure  (7)  Form  Drag  vs  Skin  Fitcdon  Drag  at  Sea  Level,  VQ=3000  fcm^ 
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Skin  Friction  Drag  (kg-m/ft^ 


Figure  (12)  Skii  Friction  Drag  V8  Width  at  See  Level,  V0»3000  kmft 
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projectile  radius  ranging  from  two  to  six  inches.  The  majority  of  the  inputs  and  the  plots 
are  in  English  units,  but  the  model  can  display  the  plots  in  SI  units  by  minor 
modifications  to  the  variables  in  the  plot  statements. 

If  the  only  consideration  is  to  maximize  the  range  and  impact  energy  of  the 
projectile,  then  minimizing  the  total  drag  will  optimize  the  projectile's  performance.  The 
form  drag  is  two  orders  of  magnitude  larger  than  the  skin  friction  drag  in  Figure  (7)  so 
the  minimum  total  drag  is  along  the  bottom  of  the  plot. 

Examining  Figure  (8)  shows  that  the  minimum  form  drag  is  the  shaved  side  ogive 
configuration  with  a  four-inch  diameter.  Refering  to  Figure  (9)  the  four  inch  wide  shaved 
side  ogive  projectile  is  forty  one  inches  long.  That  is  not  a  very  practical  length.  Figure 
(10)  or  Figure  (1 1)  shows  that  the  projectile  would  weigh  one  hundred  and  seventy-eight 
povmds. 

The  five  inch  diameter  shaved  side  ogive  projectile  has  only  slightly  higher  form 
drag  as  shown  in  Figure  (8)  and  substantially  lower  skin  friction  dr^  shown  in  Figure 
(12).  The  length  of  the  five  inch  shaved  side  ogive  projectile  in  Figure  (1 1)  is  twenty 
four  inches.  Twenty  four  inches  is  a  far  more  practical  length  for  many  reasons.  A 
twenty  four  inch  projectile  is  easier  to  store,  will  have  less  longitudinal  stresses  in  flight, 
and  the  e3q)losive  material  is  easier  to  detonate  when  it  is  packed  together  rather  than 
stretched  out.  Refering  to  Figure  (10)  or  Figure  (1 1)  the  five  inch  shaved  side  projectile 
would  weigh  one  himdred  and  forty-four  pounds. 

Using  the  same  figures  as  the  previous  two  cases,  the  six  inch  diameter  shaved 
side  ogive  projectile  is  shorter  still  at  only  fifteen  inches  in  length,  but  the  form  drag  is 
beginning  to  get  appreciably  larger.  The  skin  friction  drag  for  the  six  inch  shave  side 
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ogive  projectile  is  actually  the  minimum,  but  the  form  drag  increase  is  substantially  larger 
than  the  skin  friction  drag  decrease.  The  six  inch  shaved  wide  projectile  would  weigh 
one  hundred  and  forty-seven  pounds.  The  increased  drag  without  increased  weight 
would  mean  a  decrease  in  performance. 

Many  of  the  assumptions  made  in  defining  the  inputs  to  run  the  model  were 
selected  with  the  idea  of  replacing  a  conventional  five  inch  gun  with  a  rail  gun,  so  it 
should  be  no  surprise  that  a  modified  five  inch  projectile  is  the  apparent  best  choice.  As 
other  requirements  are  specified  the  plots  can  be  used  to  make  similar  conQ)arisons 
between  the  parameters  to  best  meet  the  desired  performance.  There  are  additional  plots 
produced  by  the  configuration  con:5)arison  program  that  weren't  needed  for  the  above 
analysis,  but  nmy  be  useful  for  other  comparisons.  Many  more  parameters  are  calculated 
in  the  model  than  are  plotted.  A  user  can  change  or  add  plot  connnands  to  get  the  plots 
necessary  to  make  the  comparisons  they  desire  for  the  requirements  they  are  given. 
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in.  STABILITY 


Projectiles  fired  over  any  appreciable  distance  require  some  form  of  stability  to 
keep  from  being  tumbled  by  minor  variations  in  wind  speed  i.e.,  wind  shear.  Projectiles 
fired  from  conventional  guns  are  ^un  in  the  barrel  so  they  have  angular  momentum  to 
provide  stability.  The  proposed  rail  gun  barrel,  being  rectangular,  will  not  impart  this 
angular  momentum  to  its  projectile.  Fins  are  the  most  practical  means  of  providing  a 
stabilizing  force. 

A  stability  program,  also  written  in  MATLAB  and  provided  in  Appendix  (C), 
models  the  stability  of  a  projectile  with  fins.  The  shaved  side  projectile  configuration 
appears  better  able  to  provide  long  range  and  high  impact  velocity  than  the  wedge  insert 
configuration  because  of  the  shaved  side's  lower  drag.  Therefore,  the  shaved  side 
configuration  is  modeled  in  the  stabihty  program.  The  equations  for  the  ogive  nose  are 
used  m  the  model,  but  if  a  spherical  nose  is  desired  it  is  simply  an  ogive  with  a  fineness 
ratio  of  one  half. 

A.  FINS 

A  variety  of  ideas  were  initially  considered  for  attaching  the  stabULzing  fins  to  the 
projectile.  Most  of  the  ideas  considered  for  attaching  fins  were  concerned  with  whether 
or  not  there  was  sufficient  room  in  the  barrel  to  accommodate  fins  large  enough  to 
provide  stability.  Attaching  large  fins  requires  that  they  be  at  least  partially  retracted 
while  in  the  barrel  of  the  rail  gun  and  that  they  deploy  after  the  projectile  has  left  the 
barrel.  Folding  fins  add  con:q)lexity  and,  as  such,  cost  to  the  over  all  system. 

Selecting  short  span  fins  with  long  chord,  as  shown  in  Figure  (13)  on  the  next 
page,  it  wdU  be  shown  will  provide  the  stability  without  adding  the  complexity  and  cost  of 
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Figure  (13)  Short  Span  Long  Chord  Fin 
Planform  View 


Front  "V^ew 

't*-  b  -« 


a  folding  fin  configuration.  A  cruciform  configuration  with  the  fins  extending  to  the 
comers  of  the  rectangular  barrel  is  probably  the  simplest  configuration  and  is  what  was 
shown  in  Figure  (1).  Symmetric  biconvex  fins  are  used  in  the  stability  program  because, 
as  per  Reference  (1),  they  are  easy  to  manufacture  and  offer  the  best  stractural  strength. 

The  dimensions  for  the  nose  and  body  are  calculated  with  the  shaved  side  ogive 
nose  equations  developed  in  the  previous  chapter.  The  dimensions  of  the  fins  are 
calculated  by  the  model  using  several  inputs  from  the  user  and  values  from  the  nose  and 
body  calculations. 

1.  Span  and  Chord 

The  span  of  the  individual  fins  is  calculated  in  equation  (3.1)  by  subtracting  the 
radius  of  the  projectile  from  the  hypotenuse  of  a  triangle  whose  sides  are  half  the 
projectile  diameter  at  the  shaved  surfeces. 

b=[(r-rf+(r-ry]*^-r  (3.1) 
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Calculating  the  span  in  this  way  keeps  the  fins  orthogonal  to  each  other  and  keeps  their 
length  correct  to  just  touch  the  barrel  wall  regardless  of  selected  projectile  radius.  Many 
of  the  aerodynamic  calculations  assume  the  fins  in  pairs  without  the  body  present  to  form 
a  clipped  delta  whose  span  is  just  twice  the  individual  fin's  span. 

The  root  chord  is  determined  as  a  user  selected  percentage  of  the  total  body  length 
calculated  for  the  selected  body  radius.  The  t^  chord  is  calculated  in  equation  (3.2)  by 
multiplying  the  span  by  the  tangent  of  the  sweep  angle  of  the  fin's  leading  edge.  A,  and 
subtracting  the  result  fi’om  the  root  chord  of  the  fin. 

Q  =  Cr  -  btanA  (3.2) 

The  selected  chord  length  has  two  effects  in  the  calculations  of  the  other  physical 
dimension.  The  longer  the  chord  length,  the  shorter  the  span  can  be  while  still  providing 
stabilization  of  the  projectile.  The  thickness  of  the  fin  is  determined  by  multiplying  the 
chord  by  a  user  inputted  thickness  to  chord  ratio,  t.  The  longer  the  chord  length,  the 
thicker  the  fin  becomes  which  increases  form  drag.  Users  are  cautioned  to  remain 
conscious  of  the  thickness  of  the  fins  as  form  drag  caused  by  uimecessarily  thick  fins  may 
severely  degrade  the  performance  of  the  projectile. 

The  sweep  angle  is  set  in  the  stability  program  to  keep  the  leading  edge  of  the  fin 
forward  of  the  Mach  cones  created  at  the  point  of  the  delta,  which  keeps  the  leading  edge 
supersonic  until  the  fi-eestream  becomes  transsonic.  If  a  user  is  sure  that  the  projectile 
will  never  slow  to  transsonic  velocities,  a  larger  sweep  angle  may  be  used  without  the 
calculations,  such  as  for  the  fin  form  drag  coefficient,  diverging.  The  sweep  angle  is  set 
in  the  fins  section  of  the  program  rather  than  in  the  input  section  to  prevent  ine3q)erienced 
users  fi'om  getting  in  trouble  by  making  a  poor  selection  of  sweep  angle. 
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2.  Surface  Areas 


The  planform  area  of  a  fin  pair  is  used  as  a  reference  area  in  many  intermediate 
calculations  relating  to  the  fins,  all  the  final  quantities  are  referenced  to  the  maximum 
body  cross  sectional  area.  The  simplest  way  to  calculate  the  planform  area  is  to  break  up 
the  fin  pair  into  a  pair  of  triangles  and  a  rectangle.  The  area  of  the  two  triangles  is  the  the 
root  chord  minus  the  tip  chord  multiplied  by  the  span.  The  rectangle  is  the  tip  chord 
multiplied  by  the  span  of  a  fin  pair.  Equation  (3.3)  combines  the  terms. 

Sf=(CrCt)b  +  2Ctb  (3.3) 

Since  the  fin  configuration  is  a  clipped  delta  wing  some  of  the  calculations  are 
based  on  what  percentage  of  the  full  delta  wing  is  present,  so  the  planform  of  the  full 
delta  wing  is  required.  The  planform  of  the  full  delta  wing  is  the  root  chord  multiplied  by 
half  the  span,  as  shown  in  equation  (3 .4) 

Sd  =  Crb  (3.4) 

3.  Volume  and  Mass 

The  thickness  of  the  fin  is  determined  from  a  thickness  to  chord  ratio,  x,  so  as  the 
chord  lessens  from  root  to  tip  so  does  the  fin  thickness.  This  makes  determining  the 
volume  of  the  fins  very  difficult.  Since  the  fins  do  not  have  a  large  sweep  angle  and  do 
have  a  very  short  span  using  the  surface  area  of  a  rectangular  fin  with  the  mean 
aerodynamic  chord  of  the  actual  fin  offers  a  very  close  approximation  to  the  surface  area 
of  the  actual  fin.  The  mean  aerodynamic  chord  for  a  clipped  delta  wing  is  calculated  with 
equation  (3.5): 

C„,a  =  (2/3)Q{[l+(Ct/Q)+(Q/Q)V[l+(Ct/Cr)]}  (3.5) 
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The  volume  can  now  be  approximated  by  determining  a  rectangular  fin  cross 
sectional  area  at  the  mean  aerodynamic  chord  and  multiplying  it  by  the  span.  To 
determine  the  fins  cross  sectional  area  the  equations  firom  Figure  (6)  can  be  used  once 
again.  The  symmetry  of  a  rectangular  biconvex  vvdng  allows  the  calculation  for  a  single 
section  to  be  multiplied  by  four  to  get  the  total  volume.  The  volume  of  the  fins  and  the 
intermediate  steps  required  to  calculate  it  are  shown  in  equations  (3.6),  (3.7),  and  (3.8). 


Rf  —  [(.5Cma)  (‘STCma)  j/T  Cma 

(3.6) 

0f  =  COS‘^(.5Craa/ Rf) 

(3.7) 

Vf  =  4b[(Rf)\7t/2  -  0f)  -  (Rf  -  .5TCma)cot0f] 

(3.8) 

The  mass  of  the  fins  can  then  be  calculated  in  equation  (3.9)  by  multiplying  the  volume 
bytwo,  for  the  two  fin  pairs,  and  the  density  of  the  component  material. 

mf=2VfPc  (3.9) 

B.  MATHEMATICAL  MODEL 

The  physical  dimensions  calculated  are  now  use  to  determine  aerodynamic 
coefficients.  The  combination  of  the  aerodyneimic  coefficients  will  be  used  to  determine 
if  the  projectile  is  stable. 

1.  Center  of  Gravity 

The  center  of  gravity  is  required  to  calculate  the  stability.  Since  the  projectile  is 
symmetric  along  its  length  the  simplest  method  of  calculating  the  center  of  gravity  is  to 
add  the  first  moments  of  all  the  components,  which  will  all  be  on  the  center  axis,  and 
divide  by  the  total  mass.  Placing  the  origin  at  the  point  of  the  nose  the  first  moment 
equations  for  the  nose,  body,  end  plate,  and  fins  become  equations  (3.10),  (3.11),  (3.12), 
and  (3.13)  respectively. 
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I„  =  .707(2fr)m„ 

(3.10) 

Ib  ~  [.5Lb+2fr][pcLb(Spcs"Secs)'*‘ItW] 

(3.11) 

Iep=  [.5(st)+Lb+2fr][pc(st)(Spcs)] 

(3.12) 

If  =  [.5C„^+(l-(%+.l))LT][2mf] 

(3.13) 

The  percent  sign  is  the  root  chord  percentage  of  body  length  and  Lt  is  the  total  projectile 
length, 

2.  Graphical  Approximations  and  Miscellaneous  Factors 

Many  aerodjmamic  coefficients  and  fectors  are  determined  empirically  and  given 
in  graphical  form.  The  four  interference  factors  from  Reference  (2)  and  the  nose/body 
lift  coefficient  from  Reference  (5)  are  presented  graphically  and  are  required  for  several 
of  the  calculations.  The  user  of  the  stability  program  would  have  to  calculate  many  of  the 
projectile's  physical  parameters  to  read  the  graphs  and  find  the  correct  inputs  to  the 
program  everytime  a  parameter,  such  as  projectile  radius  or  root  chord  percentage,  was 
changed.  Having  the  user  do  many  of  the  calculations  defeats  the  purpose  of  having  a 
computer  model.  Computers,  unfortunately,  cannot  read  graphs  so  it  became  necessary 
to  determine  equations  that  the  conq)Uter  could  use  to  calculate  the  values  of  coefficients 
and  factors.  Appendix  (D)  has  a  series  of  short  MATLAB  programs  written  to  determine 
curves  that  match  the  graphs  from  the  above  references.  The  plots  from  these  MATLAB 
programs  are  provided  in  Appendix  (D)  with  several  points,  corresponding  to  actual 
points  from  the  references,  marked  fi)r  accuracy  comparison.  Reviewing  these  plots  will 
reveal  that  the  approximations  are  seldom  more  than  a  couple  of  percent  off. 
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There  are  several  calculated  aerodynamic  quantities  that  are  used  in  the  equations. 
The  aspect  ratio  calculated  in  equation  (3 . 1 4)  is  the  square  of  the  span  divided  by  the 
planform  area  of  the  fin. 

AR  =  bVSf  (3.14) 

The  Mach  angle  p  is  the  angle  of  the  semivertex  of  the  Mach  cone  and  is  determined  in 
equation  (3.15)  by  taking  the  inverse  sine  of  the  inverse  of  the  Mach  number. 

p  =  sin'*  (1/M)  (3.15) 

The  final  quantity  is  P,  which  is  the  square  root  of  the  Mach  number  squared  minus  one, 
shown  in  equation  (3.16). 

p  =  (M^-l)*^  (3.16) 

3.  Center  of  Pressure 

The  center  of  pressure  of  the  nose/body  is  also  a  graph  determined  quantity  fi'om 
Reference  (5),  but  the  pbt  is  a  fixed  straight  line  so  long  as  the  body  length  is  at  least 
twice  the  nose  length.  Similar  to  the  nose/body  lift  coefficient,  shown  in  Appendix  (4), 
the  center  of  pressure  graph  is  divided  into  two  separate  lines  at  the  point  vsiiere  p  divided 
by  the  nose  fineness  ratio  equals  one.  An  iftthen  statement  in  the  stability  program 
determines  when  P  divided  by  the  nose  fineness  ratio  is  less  than  or  equal  to  one  and  then 
calculates  the  two  quantities  with  equations  (3.17)  and  (3.18): 

CNan  =  -.28[(p/f)-l  .9f  +  3.58  (3.17) 

XcPn  =  CG  -  [(.125LTp/f)+.15]  (3.18) 

The  center  of  pressure  is  referenced  to  the  center  of  gravity.  If  P  divided  by  the  nose 
fineness  ratio  is  greater  than  one  the  program  uses  equations  (3.19)  and  (3.20): 

CNan  =  -3.4[(ftp)-.9f  +  3.4  (3.19) 


41 


Figiire(14)  Mach  Cones  on  Fin  Planfotm 


Crossed 


Hie  efliected  area  Is 
represented  by  the 
stiipii^. 


Uncrossed 


Xcpn  =  CG  -  [(-;  125Lt£'|3)+.425]  -  (3.20) 

The  area  of  the  fin  effected  by  the  Mach  cones  created  at  the  point  of  the  delta  and 
at  the  fin  tips  is  defined  by  the  p  calculated  above.  The  effected  surface  area  will 
increase  as  the  Mach  number  decreases  and  decrease  as  the  Mach  number  increases.  The 
effected  area  has  a  lift  coefBcient  that  is  half  the  coefficient  for  the  fin  surface  in  the 
uneflfected  area  and  is  thus  an  important  intermediate  step  in  calculating  the  total  lift 
coefficient  and  center  of  pressure  location  of  the  fins.  Calculating  the  effected  sxnfece 
area  is  also  done  by  dividing  the  calculation  into  two  regions.  If  the  span  is  short  and  the 
root  chord  is  long,  as  in  Figure  (1),  the  Mach  cones  will  cross  somewhere  on  the  fin's 
surfece,  but  if  the  root  chord  gets  shorter  and  the  span  longer  or  the  Mach  number  gets 
very  large  the  Mach  cones  may  not  cross  until  after  the  fin's  trailing  edge,  as  shown  in 
Figure  (14)  above.  In  order  not  to  double  count  the  effected  surface  area  another  iftthen 
statement  is  used.  Whether  the  Mach  cones  have  crossed  before  the  trailing  edge  is 
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determined  by  adding  the  base  at  the  fin's  trailing  edge  of  the  triangles  formed  on  the 
planform  surface  by  the  Mach  cones  and  con5)aring  to  see  if  that  exceeds  the  span.  If  the 
Mach  cones  have  crossed  the  unaffected  planform  area  is  con5)uted  and  subtracted  from 
the  total  fin  planform  area  with  equation  (3.21). 

Seff  =  Sf-  {[(1.1452tan(7i/2-A-p))/(tan(7i/2-A-p)+tan(A-ir/2-p))]^  +  (3.21) 

[(1.1452tan(A-;t/2-n))/(tan(7t/2-A-p)+tan(A-7r/2-p))f  }tan(7i/2-A-p) 

If  the  Mach  cones  do  not  cross  the  calculation  of  the  effected  planform  sur&ce  area  is 

simply  the  addition  of  the  surfece  area  of  the  two  triangles  formed  by  the  Mach  cones 

shown  in  equation  (3.22). 

Seff  =  Cr^  tanp  +  tanp  (3.22) 

The  normal  force  coefficient  of  the  fins  can  now  be  determined.  If  there  were  no 
Mach  cones  the  normal  force  coefficient  would  simply  be  four  divided  by  the  quantity  p. 
The  effected  area  in  the  Mach  cones,  however,  only  has  a  normal  force  coefficient  that  is 
half  the  value  for  the  uneffected  area,  two  divided  by  the  quantity  p.  To  account  for  the 
reduced  lift  coefficient  in  the  effected  area  half  the  lift  coefficient  multiplied  by  the 
percentage  of  total  fin  planform  area  effected  is  subtracted  from  the  total  as  shown  in 
equation  (2.23). 

CN„f  =  (4/p)  -  (2Seff/pSf)  (2.23) 

The  center  of  pressure  for  a  full  delta  wing  is  simply  two  thirds  the  root  chord. 
Since  the  fins  being  used  are  basically  clipped  delta  wings  the  center  of  pressure  is 
moved  forward  and  related  to  what  percentage  of  the  full  delta  wings  area  is  actually 
present.  The  center  of  pressure  of  the  fin  is  also  effected  by  the  Mach  cones.  The  Mach 
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cone  effected  area  only  counts  as  half  its  actual  area  as  a  percentage  of  the  total  delta 
wing  area.  The  center  of  pressure  for  the  fin  is  calculated  by  equation  (3.24). 

XcPf  ~  (2/3)Cr(Sf/Sdelta  ■  •5Seff/Sdelta)  (3.24) 

The  center  of  pressure  above  is  referenced  to  the  leading  edge  at  the  root.  In  order  to 
combine  the  fin  center  of  pressure  with  the  nose/body  center  of  pressure  it  must  be 
referenced  to  the  same  point.  The  nose/body  center  of  pressure  was  easily  referenced  to 
the  center  of  gravity  since  the  tip  of  the  nose  is  always  at  the  origin  of  the  coordinate 
system.  The  position  of  the  fin’s  leading  edge,  however,  will  change  everytime  the  root 
chord  percentage  or  the  radius  of  the  projectile  is  changed.  Foreshadowing  the  drag 
calculations  of  the  next  chapter,  there  will  be  a  boattail  on  the  projectile  and  the  trailing 
edge  of  the  fin  will  be  at  the  start  of  the  boattail.  The  boattail  is  set  at  ten  percent  of  the 
body  length'so  the  leading  edge  of  the  fins  referenced  to  the  nose  of  the  projectile  will  be 
at  the  total  length  of  the  projectile  minus  the  root  chord  percentage  plus  ten  percent. 
Subtracting  the  center  of  gravity  location  from  the  fin  leading  edge  location  will  move  the 
reference  to  the  center  of  gravity  frame.  The  center  of  pressure  calculation  for  the  fins 
becomes  equation  (3.25). 

Xcpf  =  (2/3)Q(S/Sdeita  -  .5Seff/Sdeita)  -  [(1 -(%+.! ))Lt-CG]  (3.25) 

This  center  of  pressure  that  relates  the  interference  between  the  nose/body  and 
fins  and  is  found  in  graph  (15)  of  Reference  (2).  As  long  as  the  tip  chord  does  not 
become  substantially  shorter  than  the  root  chord  and  the  span  does  not  become 
substantially  larger  than  the  diameter  of  the  projectile  body  equation  (3.26)  using  a  linear 
approximation  of  the  graph  will  work. 

XcPb  =  -[  Cr(.9988P(AR)+.35)]  -  [(1-(%+.1))Lt-CG]  (3.26) 


44 


Graph  (15)  of  Reference  (2)  once  again  is  referenced  to  the  root  chord  leading  edge  so  the 
second  factor  is  used  to  rereference  the  position  to  the  center  of  gravity.  If  a 
configuration  outside  the  limits  above  is  desired  the  user  will  have  to  recalculate  the 
above  ^proximation  using  Reference  (2)  to  meet  their  requirements. 

4.  Pitching  Moment  Coefficient 

All  the  calculations  above  have  been  leading  to  the  most  important  calculation.  If 
the  pitching  moment  coefficient,  Cma,  is  negative  there  will  be  a  restoring  force  acting 
against  angle  of  attack  changes  to  keep  the  projectile  stable.  The  slope  of  the  pitching 
moment  coefficient  in  equation  (3.27)  is  the  combination  of  the  quantities,  most 
importantly  the  center  of  gravity  referenced  center  of  pressures,  calculated  above. 

CMa  =  {CNa  Xcp„  +  [CNaf  (X  cpfb  Kbw  +  X  cpf  KwB)S/2r]}/2r  (3.27) 

C.  RESULTS 

The  stability  program,  much  like  the  configuration  conq)arison  program  from  the 
previous  ch^ter,  plots  many  different  quantities  against  Mach  number.  Many  more 
combinations  are  possible,  should  the  user  desire.  The  important  plot  for  this  thesis  is  the 
plot  of  the  pitching  moment  coefficient  as  the  Mach  number  is  varied. 

The  projectile  from  Figure  (1)  with  a  long  root  chord,  fourteen  inches,  and  short 
span  about  three  quarters  of  an  inch,  and  a  five  inch  shaved  side  body  with  an  ogive  nose 
was  modeled.  This  model  produced  the  results  presented  in  Figure  (15)  on  the  next  page. 
Figure  (15)  shows  that  the  slope  of  the  pitchmg  moment  coefficient  is  indeed  negative 
through  the  Mach  region  of  interest,  so  there  will  be  a  restoring  force  against  any  angle  of 
attack  changes.  An  important  distinction  should  be  made  here.  The  fins  will  help 
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Figure  (15)  CM  alpha 


prevent  the  projectile  from  being  tumbled  by  winds  encountered  in  flight,  but  they  will 
not  prevent  the  projectile  from  being  offset  laterally  by  the  wind,  just  like  with 
conventional  projectiles.  The  next  chapter  examines  the  trajectory  of  the  projectile  and 
depending  on  the  angle  of  fire  it  is  possible  that  the  projectile  may  pass  through  the  high 
winds  of  the  jetstream.  The  speed  of  the  jetstream  winds,  normally  no  higher  than  one 
hundred  and  fifty  meters  per  second  according  to  Reference  (6),  is  small  congiared  to  the 
velocity  of  the  projectile,  normally  no  lower  than  one  thousand  meters  per  second  at 
altitude.  If  the  projectile  were  moving  at  its  slowest  and  at  ninety  degrees  to  the 
maximum  jet  stream  wind  the  angle  of  attack  would  only  be  eight  degrees.  The  angle  of 
attack  change  will  be  resisted  by  the  fins,  but  the  impact  point  of  the  projectile  may  be 
moved  substantially. 
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IV.  PERFORMANCE 

The  previous  chapters  have  provided  the  physical  dimensions  of  a  stable 
projectile  that  meets  the  requirements  of  the  rail  gun  concept.  This  chapter  reviews  the 
development  of  a  program,  written  in  the  C  computer  language  and  provided  in  Appendix 
(B),  to  model  the  performance  of  the  projectile. 

A.  PROJECTILE 

Most  of  the  dimensions  of  the  projectile  and  fins  are  calculated  with  the  equations 
from  the  configuration  comparison  and  stability  programs;  translated  into  C  for  the 
performance  program.  There  are,  however,  a  couple  of  quantities  not  previously  needed. 
The  wetted  area  of  the  fins  is  used  for  the  skin  fiiction  drag  calculation  and  the  base  area 
is  used  for  the  base  drag  calculation. 

The  wetted  area  of  the  fins  is  determined  by  using  the  same  assumptions  that  were 
used  in  determining  the  volume  of  the  fins.  The  arc  length  at  the  mean  aerodynamic 
chord  is  calculated  using  the  radius  of  curvature  of  the  fins  and  the  angle  as  determined  in 
Figure  (6).  The  wetted  surface  area  of  the  fins  is  the  product  of  the  radius  of  curvature, 
the  angle  subtended,  the  span,  and  by  four  to  account  for  the  four  such  surfaces  that  make 
up  each  fin.  Equation  (4.1)  shows  the  formula  using  the  same  notation  as  in  the  fin 
volume  calculation  of  the  previous  chapter. 

Swf=4bRfef  (4.1) 

The  base  area  of  the  projectile  is  the  area  at  the  very  end  of  the  projectile.  The 
base  area  is  reduced  from  the  projectile’s  maximum  cross  sectional  area  because  of  the 
boattail.  The  boattail  is  set  at  eight  degrees  and  as  ten  percent  of  the  total  body  length. 
Reference  (1)  shows  that  these  parameters  provide  the  best  performance  at  lower  Mach 
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numbers  which  is  where  the  projectile  will  spend  the  majority  of  its  flight.  The  radius  of 
curvature  and  the  height  of  the  rectangular  area  will  be  reduced  by  one  tenth  the  total 
length  multiplied  by  the  tangent  of  eight  degrees  as  shown  in  equations  (4.2),  (4.3),  and 
(4.4). 


rba  =  r-.014054LT 

(4.2) 

9ba  =  7i/2-sm'(h/2r) 

(4.3) 

H  ba  =  rbaC0s(sin'‘(h/2r)) 

(4.4) 

Using  equation  (2.18)  used  for  the  projectile  cross  section,  but  with  the  reduced  quantities 
from  above,  the  surface  area  of  the  base  can  now  be  determined  with  equation  (4.5). 

S  ba  =  2[(rba)  ^  Sm’(h/2r)  -  (rba)  ^COS^GbaCOteba]  - 

2hrbaCos(sin'^(h/2r))  .  (4.5) 

B.  ATMOSPHERE 

The  drag  force  on  the  projectile  is  dependent  on  the  projectile’s  environment  As 
the  projectile  moves  through  the  atmosphere  that  environment  changes,  so  the 
atmosphere,  i.e.,  the  environment,  must  be  modeled.  A  standard  atmosphere  can  be 
modeled  as  a  series  of  linear  temperature  gradients  and  isotherms  with  increasing  altitude 
as  shown  by  reference  (7).  The  temperature  gradients,  in  degrees  Kelvin  are  all  of  the 
form  given  in  equation  (4.6). 

T  =  To  +  ay/1000  ‘  (4.6) 

Where  To  is  the  temperature  at  the  bottom  of  the  gradient  in  degrees  Kelvin,  a  is  the  slope 
of  the  gradient  in  degrees  Kelvin  per  meter,  and  y  is  the  altitude  in  meters.  As  an 
example  the  To  is  two  himdred  and  eighty  eight  degrees  Kelvin  and  the  slope  is  minus  six 
and  a  half  from  zero  to  eleven  kilometers.  The  temperature  is  needed  to  determine  the 
speed  of  soimd  and  air  density  at  each  altitude. 
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The  speed  of  sound,  equation  (4.7),  at  any  altitude  is  the  square  root  of  gamma 
multiplied  by  the  gas  constant,  R,  for  air  and  the  temperature: 

a  =  (yRT)‘^  (4.7) 

The  Mach  number,  used  in  the  drag  calculations  below  one  hundred  and  five  kilometers, 
is  determined  by  dividing  the  projectile’s  velocity  by  the  local  (at  altitude)  q)eed  of 
soxmd. 

The  skin  friction  drag  of  the  projectile  at  altitude  is  directly  proportional  to  air 
density  at  that  particular  altitude.  The  air  density  for  any  altitude  is  calculated  in  equation 
(4.8)  as  a  percentage  of  sea  level  air  density  with  the  following  equation: 

Pa  =  pao(T/288)’'^^^'^**^^'‘  (4.8) 

C.  DRAG 

The  ^rformance  program  uses  six  drag  coefficients  in  determining  the  trajectory 
and  impact  velocity  of  the  projectile.  The  skin  friction  drag  and  form  drag  coefficients 
for  the  nose/body  and  fins,  and  the  base  drag  and  boattail  drag  coefficients  are  calculated 
for  the  model. 

1.  Skin  Friction  Drag 

The  generic  high  velocity  skin  friction  drag  coefficient  used  in  Chapter  n  is  not 
sufficiently  accurate  for  the  performance  calculations  because  of  the  larger  range  of 
velocities.  The  velocity  varies  from  three  thousand  meter  per  second  to  as  low  as  five 
hundred  meters  per  second.  The  projectile's  velocity  is  still  the  dominant  frictor  in  the 
numerator  of  the  Reynold's  number  equation,  equation  (2.24)  for  the  nose  body  and 
equation  (4.9)  for  the  fins. 

Ren/i)  =  pvW  |x  (2.24) 
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Ref=pvCma/M  (4.9) 

The  skin  friction  drag  coefficients  calculated  with  equation  (2.25),  using  the  appropriate 
Reynolds  number  for  the  nose/body  or  fins,  does  not  account  for  compressible  flow  at 
supersonic  speeds. 

Cf=.455  [Logio(Re)r"'^  (2.25) 

The  compressible  flow  fector  is  calculated  with  equation  (4.10). 

fcf={l  +  [(Y-l)/2]M2r""’  (4.10) 

Even  with  the  compressible  flow  fector  correction  the  skin  friction  coefficient  is  still 
close  to  .002.  Even  a  small  difference,  however,  over  a  long  distance,  can  make  a 
difference  of  kilometers  in  range,  reducing  the  accuracy  of  the  model. 

The  skin  friction  drag  coefficients  must  be  prepared  for  combination  with  the 
other  drag  coefficients  for  use  in  the  equations  of  motion  by  multiplying  them  by  the 
appropriate  wetted  area  and  dividing  by  the  reference  area,  the  projectile's  maximum 
cross  sectional  area.  These  hybrid  skin  drag  coefficients  are  the  normal  coefficients 
referenced  to  a  common  reference  area  in  order  simplify  the  computer  programming  and 
to  fecilitate  the  notation  in  this  thesis.  They  are  defined  in  equations  (4.1 1)  and  (4.12). 

Cfb  ~  Cfl)fcf(Sno/s“^Si)s)/Scss  (4.11) 

Cff’  =  CfAfS^Scss  (4.12) 

2.  Form  Drag 

The  form  drag  coefficients  are  determined  by  the  methods  from  Reference  (1). 
The  form  drag  coefficient  for  the  nose/body  is  calculated  with  the  sane  equations  as  in 
Chapter  11. 

Cpb  =  (Ap/q){  l-[(196(L„/d)2.16)/(14(M+18)(L„/d)2)]}  (2.28) 
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Ap/q  =  [.083+(.096/M^)](a/10)^-^® 

(2.29) 

a  -  tan'*(l/2f) 

(2.30) 

The  form  drag  coefficient  for  the  fins  is  equation  (4.13). 

CFf=(5.33x^Sf)/p(Sf-.5Seff) 

(4.13) 

The  biconvex  configuration  of  the  fins  accounts  for  the  5.33.  If  a  different  wing 
configuration  is  desired,  e.g.  modified  double  wedge,  the  correct  fector  must  be 
calculated  by  the  user  as  shown  in  Reference  (1)  and  changed  in  the  equations  of  motion 
as  this  factor  is  not  in  the  input  section  of  the  program.  The  terms  Sf  and  Seff  are 
determined  using  equations  (3.3)  and  (3.21)  from  Chapter  IQ. 

Sf=(Q-C,)b  +  2Qb  (3.3) 

Scff=  Sf-  {[(1.1452tan(7t/2-A-n))/(tan(7t/2-A-p)+tan(A-7i/2-p))f  +  (3.21) 

[(1.1452tan(A-7t/2-p))/(tan(7t/2-A-p)+tan(A-7c/2-n))f  }tan(7E/2-A-p) 

3.  Boattail  and  Base  Drag 

The  boattail  dr^,  Cbt,  and  base  drag,  Cb,  coefficients  are  en^irical  parameters 
presented  graphically  in  Reference  (1).  Approximations  are  calculated  in  /^endix  (D) 
with  the  other  approximation  formulas.  The  quantities  coefficients  determined  by  the 
approximations  have  to  be  adjusted  to  the  reference  area  by  dividing  by  the  maximum 
projectile  cross  sectional  area. 

4.  Drag  Factor 

The  drag  fector  used  in  the  equations  of  motion  is  a  sum  of  the  drag  coefficients 
that  have  been  referenced  to  the  projectile's  maximum  cross  sectional  area  multiplied  by 
the  dynamic  head  and  reference  area.  Equation  (4.14)  presents  the  drag  fector  with  the 
dynamic  head  shown  in  its  component  parts. 

CD  =  (Cft'+Cff +CFb+CFfi-Cbt+Cb)Scsspv^/2  (4. 14) 
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A  review  of  the  performance  program  will  reveal  that  this  equation  never  appears. 
The  factors  from  the  drag  calculations  and  equations  of  motion  are  at  times  intermingled 
in  the  performance  program  in  order  to  reduce  the  number  of  calculations  and  thus  reduce 
the  run  tune  of  the  program 
D.  EQUATIONS  OF  MOTION 

The  equations  of  motion  in  the  performance  program  were  originally  written  for 
an  assignment  in  SE  4021  involving  the  super  gun.  SE  4021  is  a  computer  simulation 
course  m  the  Physics  Department  at  the  Naval  Postgraduate  School  taught  at  the  time  by 
Dr.  WUliam  Colson.  The  class  assignment  used  a  simple  model  of  the  atmosphere,  a  set 
projectile  geometry,  and  a  set  drag  coefficient  based  on  a  "bullet"  shape.  While  this 
thesis  went  into  fer  greater  detail  actually  calculating  those  quantities,  the  actual 
equations  of  motion  are  not  changed  from  the  original  assignment. 

1.  Velocity 

The  initial  velocity  is  inputted  by  the  user  as  the  muzzle  velocity.  This  velocity  is 
decomposed  into  its  lateral  and  vertical  components  in  equations  (4.15)  and  (4.16)  by 
using  the  &ing  angle  specified  by  the  user  as  the  initial  ai^le  of  the  trajectory,  <j). 

Vx  =  vcos(j)  (4.15) 

Vy  =  vsin<j)  (4.16) 

The  performance  program  calculates  the  velocity  at  each  time  as  a  change  from 
the  previous  time  by  a  time  step.  The  time  step  is  calculated  by  dividing  the  muzzle 
velocity  by  the  acceleration  of  gravity  and  multiplying  by  a  user  selected  fector  as  in 
equation  (4.17). 

dt  =  .0001 v/g  (4-17) 
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Determining  the  time  step  in  this  fiishion  allows  the  time  step  to  be  of  reasonable  length 
for  the  selected  muzzle  velocity.  Using  a  small  factor  increases  the  number  of  time  steps, 
which  improves  the  accuracy  of  the  model,  but  slows  down  the  run  time  of  the  program. 

The  lateral  and  vertical  velocities  are  calculated  seperately  with  equations  (4.18) 
and  (4.19). 

Vx(t)  =  Vx(t-dt)  -  kv(t-dt)Vx(t-dt)dt  (4. 1 8) 

Vy(t)  =  Vy(t-dt)  -  gdt  -  kv(t-dt)vy(t-dt)dt  (4.19) 

The  vertical  velocity  is  reduced  not  only  by  the  effects  of  drag,  but  also  by  gravity.  The 
non-subscripted  velocity  is  the  total  velocity  calculated  with  equation  (4.20). 

v  =  (vx"  +  vy2)i^  (4.20) 

Ilie  fector  k  is  related  to  the  fector  CD,  calculated  in  the  previous  section,  by  equation 
(4.21). 

k  =  CD/mv^  =  (Cfl,'+Cff'+CFb+CFi+Cbt+Cb)Scssp/2m  (4.21) 

Using  k  instead  of  CD  in  the  equations  of  motion  greatly  simplifies  the  formulas. 
Dividing  CD  by  m  converts  the  drag  force  into  an  acceleration.  Multiplying  an 
acceleration  by  a  time  step  provides  a  velocity.  Dividing  CD  by  the  velocity  squared  is  a 
contrivance  to  avoid  having  to  track  the  angle  of  the  trajectory  at  each  time  step.  In 
equations  (4.18)  and  (4.19)  the  v^,  which  came  from  the  dynamic  head  calculation,  is 
replaced  by  v(t)vx(t-dt)  and  v(t)vy(t-dt),  as  £q)propriate,  which  have  the  angle  information 
imbedded  in  the  subscripted  velocity  as  developed  for  the  lateral  case  in  equation  (4.22). 

Vx  =  vcoscj)  .’.  v^cos(j)  =  v(vcos(|))  =  Wx  (4.22) 
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Altitude  (km) 


2.  Position 


The  lateral  and  vertical  positions  are  once  ^ain  calculated  separately  with 
equations  (4.23)  and  (4.24)  by  adding  the  change  in  position  since  the  previous  time  step 
to  the  previous  time  step's  position. 

x(t)  =  x(t-dt)  +  Vxdt  (4.23) 

y(t)  =  y(t-dt)  +  Vydt  (4.24) 

E.  RESULTS 


The  projectile  has  a  maximum  range  of  one  hundred  and  eighty  two  kilometers,  as 
seen  in  Figure  (16)  above,  when  fired  with  a  fifty  seven  degree  launch  angle.  That  range 
is  approximately  eight  times  the  maximum  range,  twenty  two  kilometers,  of  a 
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conventional  five  inch  gun.  Figure  (17)  below  shows  that  the  projectile  only  takes  four 
minutes  and  eight  seconds  to  reach  a  target  at  its  maximum  range  and  eight  seconds  to 
reach  a  target  at  the  maximum  range  of  a  conventional  five  inch  gun. 


Figure  (17)  Impucl  Time 


The  kinetic  energy  of  the  projectile  at  impact  is  also  substantially  larger  than  a 
conventional  five  inch  projectile.  Figure  (18)  on  the  next  page  shows  the  energy  of  the 
projectile  at  impact  for  ranges  from  zero  to  maximum.  There  are  two  energies  for  almost 
all  the  ranges  because  the  launch  velocity  is  held  constant  and  the  launch  angles  are 
varied  between  zero  and  eighty  nine  and  a  half  degrees.  An  interesting  feature  of  Figure 
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(18)  is  that  for  ranges  greater  than  sixteen  kilometers  the  projectile  has  greater  kinetic 
energy  if  shot  at  the  high  angle  rather  than  the  direct  angle. 


Two  mechanisms  are  at  work  to  improve  the  performance  of  the  high  angle  shots. 
The  first  mechanism  is  that  the  high  angle  shots  reach  the  higher  atmosphere,  which  has 
lower  density  and  thus  less  drag,  quicker  so  they  do  not  slow  down  as  rapidly.  The 
second  mechanism  is  that  the  higher  angle  of  fire  shots  also  descend  at  a  higher  angle  so 
gravity  helps  recoop  some  of  the  drag  loss,  i.e.  the  vertical  component  of  velocity  is  a 
larger  part  of  the  total  velocity. 
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Some  of  the  points  of  interest  on  Figure  (18)  are  the  maximum  energy,  the  energy 
CTOssover  point  between  direct  fire  and  high  angle  fire,  the  minimum  energy,  and  the 
energy  at  maYimutn  range.  At  the  muzzle  of  the  rail  gxm  the  projectile  has  a  maximum 
kinetic  energy  of  three  hundred  and  seven  megajoules,  this  is  equivalent  to  the  energy 
released  by  detonating  two  hundred  and  forty  four  pounds  of  TNT.  The  crossover  point 
where  the  kinetic  energy  of  a  projectile  fired  at  a  shallow  angle  (direct  fire)  and  a  high 
angle  occurs  at  about  sisteen  kilometers.  The  kinetic  energy  at  sixteen  kilometers  is 
twenty  seven  megajoules,  this  is  equivalent  to  twenty  one  pounds  of  TNT.  Referring  to 
Figure  (17),  the  projectile  takes  either  eleven  or  three  himdred  and  twenty  five  seconds  to 
arrive.  The  minitnum  kinetic  energy  of  six  point  eight  megajoules,  equivalent  to  only 
five  and  a  half  pounds  of  TNT,  occurs  at  fidly  kilometers  with  a  launch  angle  of  twenty 
one  and  a  half  degrees.  The  projectile  reaches  fifty  kilometers  in  one  minute.  As 
mentioned  above  the  projectile  reaches  its  maximum  range  of  one  hundred  and  eighty 
two  kilometers  in  four  minutes  and  eight  seconds.  The  projectile  has  eighteen 
magajoules  of  kinetic  ei^rgy  when  it  impacts  at  its  maximum  range,  this  energy  is 
equivalent  to  fifteen  pounds,  of  TNT. 

The  projectile  often  has  the  kinetic  energy  equivalent  to  an  amount  of  TNT 
substantially  larger  than  the  ejqplosive  payload  of  the  projectile.  The  destructive  power, 
however,  may  be  substantially  less  than  the  TNT  because  there  is  no  mechanism  to 
efiBciently  transfer  the  projectile's  kinetic  energy  to  the  target. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

An  aerodynamically  stable  projectile  can  be  made  for  a  rectangular  barrel  Figure 
(1)  sbows  the  basic  configuration  and  dimensions  of  the  projectile  developed  and  then 
used  as  an  example  in  this  thesis.  The  projectile  performs  very  well,  but  the  performance 
model  is  for  the  best  e3q)ected  muzzle  velocity  and  with  a  mass  that  may  be  too  large  for 
practical  shipboard  deployment.  The  real  benefit  of  this  thesis  is  the  three  computer 
programs  that  will  allow  future  users  to  quickly  analyze  different  projectile 
configurations  to  meet  evolving  requirements. 

A  great  deal  can  be  done  to  improve  the  capabilities  of  the  programs  developed  in 
this  thesis.  The  development  and  incorporation  of  the  equations  for  the  subsonic  regime 
into  the  models,  with  a  transsonic  approximation  to  link  the  subsonic  and  supersonic 
regimes,  will  allow  a  greater  variety  of  materials  to  be  considered.  Different  fin  types, 
i.e.  double  wedge  or  modified  double  wedge,  can  also  be  incorporated  into  the  models  to 
give  USCTS  a  choice. 

A  proper  heat  transfer  study  of  the  projectile  should  be  done.  The  information 
fi'om  a  heat  transfer  study  can  be  used  to  determine  both  the  fineness  ratio  at  which  the 
ogive  nose  would  melt  and  what  materials  can  be  used  in  the  construction  of  the 
projectile.  The  form  drag  of  the  nose  is  dependent  on  the  fineness  ratio.  The 
performance  of  the  projectile  can  be  improved  by  making  the  nose  pointier  to  reduce  the 
drag,  but  if  the  nose  is  too  pointy  it  may  melt.  Different  materials  may  be  used  for 
different  portions  of  the  projectile  if  the  tenq)erature  of  a  desired  conqwnent  is  known  not 
to  exceed  the  melting  point  of  the  con^nents  constituent  material. 
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Wind  tunnel  testing  could  verity  the  drag  equations  used  or  provide  information 
to  improve  their  accuracy.  Hypersonic  wind  tunnels  are  very  expensive  to  operate,  but 
the  proposed  projectile  flies  most  of  its  trajectory  at  velocities  easily  achieved  in  standard 
wind  tunnels.  A  study  of  the  true  effects  of  the  flat  contact  surfaces  would  determine  if 
the  assumptions  made  in  calculating  the  drag  of  the  projectile's  body  sufficiently 
accounted  for  the  airflow  separation  and  turbulence  that  the  flat  contact  surfaces  will 
probably  cause. 

The  most  difficult  study  required  may  well  be  determining  what  explosive 
compound  to  use  in  the  projectile  and  how  it  can  be  fused.  The  ejqilosive  conipound 
used  will  have  to  be  able  to  undergo  the  shock  of  the  rapid  acceleration  experienced  in 
the  barrel  of  the  rail  gun  and  the  high  temperatures  the  projectile  may  reach  without 
detonating.  Electrical  foses  may  be  damaged,  or  worse,  detonated  by  the  very  high 
currents  encountered  in  the  barrel.  An  impact  fuse  in  the  base  may  work,  but  a  thorough 
study  should  be  done  to  insure  the  projectile  does  not  detonate  in  the  barrel  of  the  rail 
gun. 

Experiments  with  rail  guns  are  on  going.  This  thesis  provides  tools  to  help  in  the 
design  of  the  total  weapon  system. 
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APPENDIX  (A)  MATLAB  MODEL  FOR  NOSE/BODY 
CONFIGURATION  COMPARISON 


%  Inputs 


stin=.5;  %  skin  thickness  in  inches 

rhos=19200.0;  %  density  of  metal  in  kg/m^3  (tungsten) 

rhoe=1520-0;  %  density  of  explosive  in  kg/m-'B  (04) 

v0=3000;  %  initial  shell  velocity  in  m/s 

cwin=2;  %  width  of  contact  area  in  inches 

ew=5.5;  %  explosive  weight  in  kg 

f=l;  %  fineness  ratio,  (nose  length) /diameter 

%  Constants 

rho=1.03;  %  ai'r  density  at  sea  level  in  kg/m''3 

a=346;  %  speed  of  sound  in  air  at  seal  level/standard  day  in  m/s 

%  Calculated  Initial  Values 

M0=v0/a; 

q=.5*rho*  (v0'‘2) ;  %  kg/m-s^2 

stm=stin*. 025416;  %  skin  thickness  in  meters 

cwm=cwin*. 025416;  %  width  of  contact  area  in  meters 

%  Drag  Factors 

Cf=.002;  %  skin  friction  drag  coefficient  from  AA3710  notes  pg  6-4 
sigmasph=45;  %  semivertex  angle  in  degree  for  sphere 

sigmaogi=57.2958*atan(l/ (2*f) ) ;  %  semivertex  angle  in  degrees  for  ogive 
presssph=  ( .  083+  ( .  096/  (M0^2) ) )  *  ( (sigmasph/10)  -'1 . 69)  ; 

%  from  AA3710  notes  pg  12-2 
pressogi=  ( .  083+  ( . 096/  (M0''2) ) )  *  ( (sigmaogi/10)  ^1.69)  ; 

CDFsph=presssph* (1- ( ( (392* ( . 5^2) ) -32) / (26* (MO+18) * ( .5^2) ) ) ) ; 
CDFogi=pressogi* (1- ( ( (392* (f-2) ) -32) / (26* (MO+18) * (f^2) ) ) ) ; 

%  Spherical  with  Wedge  Insert 

srinw=Linspace(2,6,30) ;  %  radius  in  inches 
sdinw=(2.*srinw)+cwin;  %  Projectile  width  in  inches 
srmw=srinw*. 0254166667;  %  radius  in  meters 
snvw=(  ( (2/3)  *pi)  .  *  (srmw. -^3) )  +  (  ( (pi/2)  *cwm)  .  *  (srmw. ''2) ) ; 

%  Nose  volume  in  m''3 

snww=rhos.*snvw;  %  Weight  of  nose  in  kilograms 

spcsw=(pi.*  (srmw.''2) )  +  ( (2*cwm)  .*srmw)  ; 

%  Projectile  cross  section  in  meters ''2 
spcsiw— spcsw ♦ / • 00 0 64 6 ;  %  Projectile  cross  section  in  in  2 
secsw=(pi.*(  (srmw-stm)  . ''2) )  +  ( (2*cwm)  .*(srmw-stm) ) ; 

%  Explosives  cross  section  in  m^2 
selw=ew./(rhoe.*secsw) ;  %  Length  of  explosive  in  meters 
spwkw=(rhos-* ( ( (spcsw-secsw) . *selw) + (stm. *spcsw) ) )+ew+snww; 

%  Projectile  weight  in  kilograms 
spwlw=spwkw.*2.2;  %  Projectile  weight  in  Ibm 
stlmw=srmw+selw+stm;  %  Total  length  of  projectile  in  meters 
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stliw=stlmw/- 0254166667;  %  Total  length  in  inches 

ssamw=  (2 .  *pi. *  {srmw.  ^2)  )  +  {pi.*srmw. *cwm)  +(2 .  *pi.  *srmw.  *  (selw+stm)  )  .  .  . 

+  (2 . * (selw+stm) . *cwm)  ;  %  Projectile  Surface  Area  in  m^2 

ssaiw=ssamw. /. 000646;  %  Projectile  Surface  Area  in  in^2 
ssfdmw=Cf *ssamw*q;  %  skin  friction  drag 
sfdmw=CDFsph*spcsw*q;  %  form  drag 

%  Spherical  with  Shaved  Sides 

srins=Linspace (2 . 0, 6, 30) ;  %  radius  in  inches 

sdins= (2 . *srins) ;  %  Projectile  width  in  inches 

srms=srins*. 0254166667;  %  radius  in  meters 

snvs=(  (  (2/3)  *pi)  .  *  (srms.  ^3) )  -  (  (cwm.  *  (srms,  "^2)  *pi.*.5)  *  . 

(cwm.  *  (srms .  ^"2)  .  *acos  (cwm.  /  (4.  *srms) ) )  -  (4 .  *  (srms .  ^3)  .  *  (sqrt  (1-. , . 

(  ( (cwm./(4.*srms)  )  .^^2)  )  ) )  )  +  (4.*  (srms,'^3) )~(  (srms . *cwm. /24)  . 

{ (cwm. / (2 . *sin (acos (cwm. / (16. *srms) ) ) ) ) + (cwm. / (2 . *sin (acos (cwm. / . . . 
(8.*srms) ) ) ) )+(3.*cwm./ (2.*sin(acos { (3.*cwm) ./(16.*srms) ))))+... 
(cwm. / (2. *sin{acos (cwm. / (4 . *srms) ))))))  +  ( (cwm./ (96. *srms)  ).*... 

( ( (cwm. ^3) . / (128 . *sin (acos (cwm. / (16. * srms) ) ) ) ) + ( (cwm. ^3) . / (32 . * . . . 

sin (acos (cwm. / (8 . *srms) ))))  +  ( (27.*  (cwm. ^3) ) ./ (128 .*sin (acos ( (3.*cwm)  . . . 
. / (16. *srms) ))))+( (cwm. ^3) . / (8 . *sin (acos (cwm. / (4 . *srms) ))))))); 

%  Nose  volume  in  m^3 

snws=rhos.*snvs;  %  Weight  of  nose  in  kilograms 
spcss=  (pi.  *  (srms  -""2)  )  -  (2  .  *  (  (  (srms,  ^^2)  .  *  (  (pi/2)  . . 

acos  ( .  5*cwm.  /srms)  )  )  -  (  (  (srms-  (srms .  *  (1- 
sin (acos (,5.*cwm./srms) ))))... 

. ^2) . *cot (acos { , 5*cwm. /srms) ) )  )  )  ; 

%  Projectile  cross  section  in  meters^2 
spcsis=spcss ,/. 000646;  %  Projectile  cross  section  in  in^2 
secss= (pi. * ( (srms-stm) . ^2) ) - (2 . * ( ( ( (srms-stm)  . ^2) . * ( (pi/2) -. . . 
acos  ( .  5*cwm.  /  (srms-stm)  )))-(((  (srms-stm)  -  (  (srms-stm)  .  * .  .  . 

(1-sin  (acos  ( .5 .  *cwm.  /  (srms-stm)  ))))).  ^^2)  .  * ..  . 

cot (acos ( .5*cwm. / (srms-stm)  ))))) ;  %  Explosives  cross  section  in  m^2 

sels=ew. / (rhoe. *secss) ;  %  Length  of  explosive  in  meters 
spwks= (rhos . * ( ( (spcss-secss) . *sels) + (stm. *spcss) ) ) +ew+snws; 

%  Projectile  weight  in  kilograms 
spwls=spwks.*2.2;  %  Projectile  weight  in  Ibm 
stlms=srms+sels+stm;  %  Total  length  of  projectile  in  meters 
stlis=stlms/. 0254166667;  %  Total  length  in  inches 

ssams= (2 .*pi.* (srms. ^2) ) - (2. *pi. * ( ( (- (srms .^2) . /2) .*log( (.5. *cwm) + _ 

sqrt ( (srms. ^2) +(.25.* (cwm^2) ) ) ) ) + ( . 25. *cwm. *sqrt ( (srms. ^2) + ( .25. * . . . 
(cwm^2) ) ) )  +  ( .5. * (srms. ^2) . *log (srms) ) ) )  +  (pi.* ( (cwm. *.5) .^2) )  +  . . . 

(4 . *srms.* (sels+stm) . *asin ( . 5. *cwm. /srms) )  +  (2 . *cwm. * (sels+stm)  )  ; 

%  Projectile  Surface  Area  in  m'^2 
ssais=ssams. /. 000646;  %  Projectile  Surface  Area  in  in^2 
ssfdms=Cf *ssams*q;  %  skin  friction  drag 
sfdms“CDFsph*spcss*q;  %  form  drag 


%  Ogival  with  Wedge  Insert 

orinw=Linspace (2, 6^ 30) ;  %  radius  in  inches 
odinw=(2.*orinw) +cwin;  %  Projectile  width  in  inches 
ormw=orinw. *.0254166667;  %  radius  in  meters 
onlw=2 . *ormw. *f ;  %  nose  length  in  meters 
oRmw=( (onlw. ^2) + (ormw. ^2) ) ./(2.*ormw) ; 
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%  Radius  of  curvature  of  nose  in  meters 
thetaw=acos (onlw. /oRmw) ; 

onvw={pi.*(  (onlw.*  ( (2.*  {oRmw.''2) )- (2.*oBmw.*ormw)  +  (ormw.''2) ))-... 

( (onlw.  ''3)  .  /3)  -  (onlw.  *  (oRmw-ormw)  .*sqrt  { (oRmw.^2)  -  (onlw.  ''2) ) )  *  •  • 

-  ( (oRmw-ormw)  .*  (oRmw.^2)  .  *asin(onlw.  /oRmw) ) )  )  +  (cwm.*  ( ( (oRmw.  ''2)  .  . . 
.*  (  (pi/2) -thetaw) )-( ( (oRmw-ormw)  .''2)  . *cot (thetaw) ) ) ) ; 

%  Nose  volume  in  m^3 

onww=rhos . *onvw;  %  Weight  of  nose  in  kilograms 
opcsw=(pi.*{ormw.^2) )+(2.*cwm.*ormw) ; 

%  Projectile  cross  section  in  meters^2 
opcsiw=opcsw. /. 000646;  %  Projectile  cross  section  in  in  2 
oecsw=  (pi.  *  ( (ormw-stm)  . ^'2) )  +  {2 .  *cwm. *  (ormw-stm) ) ; 

%  Explosives  cross  section  in  m''2 
oelw=ew./(rhoe.*oecsw) ;  %  Length  of  explosive  in  meters 
opwkw= (rhos . * ( ( (opcsw-oecsw) . *oelw) + (stm. *opcsw) ) ) +ew+onww; 

%  Projectile  weight  in  kilograms 
opwlw=opwkw.*2.2;  %  Projectile  weight  in  Ibm 

otlmw=(f .*2.*ormw) +oelw+stm;  %  Total  length  of  projectile  in  meters 
otliw=otlmw/. 0254166667;  %  Total  length  in  inches 
osamw=(2.*pi.* ( (2.*f.*ormw.*oRmw)-(oRmw.* (oRmw-ormw) .*asin. . . 

(2.*f .*ormw./oRmw) ) ) )+(2.*cwm.*oRmw.* ( (pi/2) -thetaw) )+. . . 

(2 . *pi. *ormw. * (oelw+stm) )  +  (2 . *cwm. * (oelw+stm) ) ; 

%  Projectile  Surface  Area  in  m''2 
osaiw=osamw./. 000646;  %  Projectile  Surface  Area  in  in ''2 
osfdmw=Cf*osamw*q;  %  skin  friction  drag 
ofdmw=CDFogi*opcsw*q;  %  form  drag 


%  Ogival  with  Shaved  Sides 

orins=Linspace(2.0,6,30) ;  %  radius  in  inches 
odins=(2.*orins) ;  %  Projectile  width  in  inches 
orms=orins*. 0254166667;  %  radius  in  meters 
onls=2.*orms.*f;  %  nose  length  in  meters 
oRms=(  (onls.''2)  +  (orms.''2) )  ./  (2.*orms) ; 

%  Radius  of  curvature  of  nose  in  meters 
thetas=acos (onls. /oRms) ; 

onvs=(pi.*  ( (onls.*  ( (2 . *  (oRms. ^2) )  -  (2 . *oRms .  *orms)  +  (orms. ^^2) ))-... 

( (onls.''3)  ./3)-(onls.*(oRms-orms)  .*sqrt  ( (oRms  .''2)  -  (onls.  ^2) ) )  •  •  • 

-  ( (oRms-orms)  .*  (oRms. ^2)  .*asin  (onls . /oRms) ) ) )  -  ( (cwm.*  (oRms.  ^"2)  * •  •  • 
pi.*.5)-(cwm.*(oRms.''2)  .*acos  ( (f  .*cwm)  .  /  (2 .  *oRms) ) )  -  ( (2 . /f )  .*... 
(oRms. ^3) .*(sqrt(l-(((f.*cwm) . / (2 . *oRms) ) . ^2) ) ) ) + ( (2 . * (oRms . ^ . . . 

3) ) ./f ) - ( (f .*oRms. *cwm. /12) . * ( (cwm. / (2 . *sin (acos ( (f . *cwm)  . / (8 . * . . . 
oRms) ) ) ) )+(cwm./ {2.*sin(acos ( (f .*cwm) ./ (4.*oRms) ) ) ) )+(3.*cwm./. . . 

(2 .*sin(acos ( (3.*f . *cwm) ./ (8. *oRms) ) ) ) )  +  (cwm.  / (2.*sin(acos ( (f .*  — 
cwm)  ./(2.*oRms)))))  ))  +  {({cwm.*(f."3))  ./(12.*oRms))  .*(((cvnn.''3)  ./... 
(128 .*sin (acos ( (f .*cwm) . / (8 . *oRms) ))))+( (cwm. ^3) ./ (32.*sin(acos — 
((f.*cwm)  ./(4.*oRms)))) )  +  ((27.*(cwra.-3) ) . / (128 . *sin (acos { (3 . *f . * . . . 

cwm) . / (8.*oRms) ))))  +  ( (cwm. ^3) . / (8 .*sin (acos ( (f .*cwm)  ./ (2.*oRms) ))))))) 
%  Nose  volume  in  m''3 

onws=rhos . *onvs ;  %  Weight  of  nose  in  kilograms 

opcss=(pi.*(orms.''2)  )-(2.*(  ((orms.''2)  .*(  (pi/2)-acos(.5*cwm./orms) ) )-. . 

( ( (orms-  (orms.  *  (1-sin  (acos  ( .  5 .  *cwm.  /  orms) )}))  .''2)  .  *cot  (acos  ( .  5 .  * . . . 
cwm./orms) ) ) ) ) ; 

%  Projectile  cross  section  in  meters''2 
opcsis=opcss./. 000646;  %  Projectile  cross  section  in  in''2 


63 


oecss=  (pi.^  ( (orms-stm)  .  ^2)  )  -  (2.  *  (  {  (  (orms-stm)  .  ^2)  .  *  ( {pi/2)  - .  •  • 
acos  ( -5*cwm.  /  (orms-stm)  )))-{(( (orms-stm)  -  ( (orms-stm)  .  * .  . , 
(l-sin(acos{.5.*cwm./  (orms-stm) ))))).  ^^2)  • 

cot (acos ( .5*cwm. / (orms-stm) ))))) ;  %  Explosives  cross  section  in  m^2 

oels=ew. / (rhoe, *oec5s) ;  %  Length  of  explosive  in  meters 
opwks= (rhos . * ( ( (opcss-oecss)  . *oels)  +  (stm, *opcss) ) ) tew+onws; 

%  Projectile  weight  in  kilograms 
opwls=opwks • *2 .2;  %  Projectile  weight  in  Ibm 

otlms-(f .*2.*orms) +oels+stm;  %  Total  length  of  projectile  in  meters 

otlis=otlms/. 0254166667;  %  Total  length  in  inches 

osams- (2,*pi.*  (  (2 .  *f .  *orms .  *oRras)  7  (oRms .  *  (oRms-orms)  ,,*asin.  •  , 

(2 .  *f .  *orms ♦  /oRms)  )  )  )  -  (oRms .  *  (  (pi.  *cwm)  -  (2 .  *cwm.  *acos  (  (f .  *cwm)  .  /  . .  . 
{2.*oRms) ) )-(  {4.*oRms.*sqrt(l-(  (  (f.^cwm)  ./{2.*oRms)  )  . ^2)  )  )  . /f ) +  . . . 

(4 .  *oRms .  ft) )  }  +  (2 .  *  ( {  (oRms .  ^^2)  .  *  ( (pi.  /2)  -acos  (cwm.  /oRms)  ))-... 

( (  (oRms-  (cwm/2)  )  .  ^2)  .  *cot  (acos  (cwm.  /oRms) ) ) ) )  +  (4 .  *orms.  *  {oels+.  - . 
stm)  .  *asin  ( .  5 .  *cwm.  /orms)  )  +  (2  .*cwm.  *  (oels+stm) )  ; 

%  Projectile  Surface  Area  in  m'^2 
osais=osams./. 000646;  %  Projectile  Surface  Area  in  in^2 
osfdms=Cf*osams*q;  %  skin  friction  drag 
ofdms=CDFogi*opcss*q;  %  form  drag 

figure  (1) 

plot  ( sdinw,  spcsiw,  ’  k-  \  sdins ,  spcsis ,  *  k .  * ) ,  grid 
title (* Cross  Section  vs  Width*) 

xlabel (’ Projectile  Width  (in) *), ylabel (* Projectile  Cross  Section 
(in-2) M 

text (6. 1, 110, *  Solid  =  w/  Insert*) 
text (6.1, 105, *  Dotted  =  w/  Shaved  Sides’) 
figure (2) 

plot  ( sdinw,  spwlw,  *  k-  * ,  sdins ,  spwls ,  *  k .  * ,  odinw,  opwlw,  *  k .  - .  * ,  odins ,  .  .  . 

opwls ,  *  k —  ’ )  ,  grid 
title { ’Weight  vs  Width’) 

xlabel (’ Projectile  Width  (in)’),  ylabel (’ Projectile  Weight  (Ibm)’) 

text (6. 1, 654, ’ Solid  =  Spherical  w/  Insert’) 

text (6. 1, 614,  ’ Dotted  =  Spherical  w/  Shaved  Sides’) 

text (6. 1, 574, ’Dash/Dot  =  Ogival  w/  Insert’) 

text (6. 1, 534, ’Dashed  =  Ogival  w/  Shaved  Sides’) 

figure (3) 

plot  ( sdinw,  s tliw,  ’  k-  ’ ,  sdins ,  stlis ,  ’  k .  ’ ,  odinw,  otliw,  ’  k  -  - .  ’ ,  odins ,  .  . . 

otlis, *  k —  * ) , grid 
title ( ’Length  vs  Width’) 

xlabel  (’ Projectile  Width  (in^)’),  ylabel  ( ’Total  Length  (in)’) 

text (6. 1, 41, ’Solid  =  Spherical  w/  Insert*) 

text (6.1,38.5, ’Dotted  =  Spherical  w/  Shaved  Sides’) 

text ( 6. 1, 36, ’ Dash/Dot  =  Ogival  w/  Insert’) 

text(6.1, 33.5, ’Dashed  =  Ogival  w/  Shaved  Sides’) 

figure  (4) 

plot (spwlw, stliw, ’ k- ’ , spwls, stlis, ’ k, ’ , opwlw, otliw, ’ k. -. * , opwls, . . . 

otlis, ’ k — ’ ) , grid 
title ( ’Length  vs  Weight’) 

xlabel (’ Projectile  Weight  (Ibm)’),  ylabel ( ’Total  Length  (in)’) 

text  (205, 41,  *  Solid  ==  Spherical  w/  Insert’) 

text (205,  38 . 5,  ’ Dotted  =  Spherical  w/  Shaved  Sides’) 

text (205, 36, ’Dash/Dot  =  Ogival  w/  Insert’) 

text (205,  33.5, ’Dashed  =  Ogival  w/  Shaved  Sides’) 

figure (5) 

plot (spwlw, spcsiw, ’k-’ , spwls, spcsis, ’k. ’ , opwlw, opcsiw, ’k.-. ’ ,  opwls,  . . - 
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opcsis, ' k — ' ) , grid 
title ('Cross  Section  vs  Weight') 

xlabel ('Projectile  Weight  (Ibm) '), ylabel (' Projectile  Cross  Section 
(in-2) ') 

text (215,  37, 'Solid  =  Spherical  w/  Insert') 
text (215, 30, 'Dotted  =  Spherical  w/  Shaved  Sides') 
text (215,23, 'Dash/Dot  =  Ogival  w/  Insert') 
text (215, 16, 'Dashed  =  Ogival  w/  Shaved  Sides') 
figure (6) 

plot (spcsiw, ssaiw, ' k- ' , spcsis, ssais, ' k. ' , opcsiw, osaiw, ' k. - . ' , opcsis, . . . 
osais, ' k — ' ) , grid 

title ('Surface  Area  vs  Cross  Section') 
xlabel ( 'Projectile  Cross  Section  (in-2)') 
ylabel ('Projectile  Surface  Area  (in-2)') 
text (21,  480, 'Solid  =  Spherical  w/  Insert') 
text (21, 460, 'Dotted  =  Spherical  w/  Shaved  Sides') 
text (21, 440, 'Dash/Dot  =  Ogival  w/  Insert') 
text (21, 420, 'Dashed  =  Ogival  w/  Shaved  Sides') 
figure (7) 

plot(ssfdmw,sfdmw./le5,'k- 

'  ,ssfdms,sfdms./le5,  'k. ' ,  osfdmw,  ofdmw. /Ie5, . . . 

'k.-. ' , osfdms, ofdius . /Ie5,  'k — ' )  ,grid 
title ('Form  Drag  vs  Skin  Friction  Drag  at  Sea  Level,  V0=3000  km/s') 
xlabeK’Skin  Friction  Drag  (kg-m/s -2) '),  ylabel  (' Form  Drag  x  100,000 
(kg-m/s-2) ' ) 

text(1750, 1.56, 'Solid  =  Spherical  w/  Insert') 
text (1750, 1.41, 'Dotted  =  Spherical  w/  Shaved  Sides') 
text (1750, 1.26, 'Dash/Dot  =  Ogival  w/  Insert') 
text(1750, 1.11, 'Dashed  =  Ogival  w/  Shaved  Sides') 
figure (8) 

plot(sdinw, ssfdmw, ' k- ' , sdins, ssfdms, 'k. ' ,odinw, osfdmw, 'k.-. ',odins, . . . 
osfdms, ' k — ' ) , grid 

title ('Skin  Friction  Drag  vs  Width  at  Sea  Level,  V0=3000  km/s') 

xlabel ( 'Projectile  Width  (in) '), ylabel (' Skin  Friction  Drag  (kg-m/s-2)') 

text (6. 1,2900, 'Solid  =  Spherical  w/  Insert') 

text (6. 1,2800, 'Dotted  =  Spherical  w/  Shaved  Sides') 

text (6. 1,2700, 'Dash/Dot  =  Ogival  w/  Insert') 

text (6.1, 2600, ' Dashed  =  Ogival  w/  Shaved  Sides') 

figure (9) 

plot (sdinw, sfdmw. /Ie5, ' k- ' , sdins, sfdms . /Ie5, ' k. ' , odinw, ofdmw . /Ie5, . . . 

' k. - . ' , odins, ofdms . /Ie5, ' k — ' ) ,  grid 
title (Form  Drag  vs  Width  at  Sea  Level,  V0=3000  km/s') 

xlabel ('Projectile  Width  (in) ' ) ,ylabel ( ' Form  Drag/100,000  (kg-m/s-2)') 

text (5, 2. 41, 'Solid  =  Spherical  w/  Insert') 

text (5, 2. 2 4, 'Dotted  =  Spherical  w/  Shaved  Sides') 

text (5, 2. 07, 'Dash/Dot  =  Ogival  w/  Insert') 

text (5, 1.90, 'Dashed  =  Ogival  w/  Shaved  Sides') 


If  a  color  printer  is  available  different  colored  lines  are  a  better  presentation  on  the  plots. 
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APPENDIX  (B)  PERFORMANCE  PROGRAM  (in  C) 

/♦Performance*/ 

#include  <stdio.h> 

#include  <inath.h> 

main(){  ^  « 

float  g,  C,  rhoo,  rho.  A,  v,  vx,  vy,  x,  y,  theta,  t,  dt,  k,  h,  xf ,  yf , pi,  r,  f ; 
float  a,M, T,  Reb, Ref,  CDFb, Cfb, CDfb, pr,  gamma,  R, mu,  SA,  pb, dc,  dee, er ; 
float  CDbt , CDb , nb , lb , Sb , r in , s tin , r hom,  r hoe ,  s tm,  nl , nv , nw, br , lb t ; 
float  ecs r el, ew, nlm, nlm2 , nlm3, tlm,  1,  sigma, ein, r2, r3,  sigmac,  ewin; 
float  cr,  ct,lamdam,mac,  tau,  tmac,rf  ,phi,bin,b,  fv,  fw,  fsa,m,  cw3; 
float  ibrh, ibrw, omega, pf , Cff, CDff, CDFf , beta, Mu, sef f , sf , et, nbt; 
float  lamda, psi, Psi, Dl, D2, D3, comf , rn, rn2, rn3, thetan, crp, f 3, cw; 
int  j; 

/*  INPUTS  */ 

rin=2.5;  /*  radius  in  inches  */ 
stin=.5;  /*  skin  thickness  in  inches  */ 

■  cwin=2;  /*  width  of  contact  area  in  inches  */ 
rhom=19200.0; 

/*  density  of  metal  in  kg/m^3  (tungsten  in  example)  */ 
rhoe=1520.0; 

/*  density  of  explosive  in  kg/m''3  (C4  in  this  example)  */ 
ew=5.5;  /*  desired  explosive  payload  in  kg  */ 
v=3000.0;  /*  muzzle  velocity  in  m/s  */ 

theta=. 986111;  /*  launch  angle  in  radians  (56.5  degrees)  */ 
tau=^.03';  /*  fin  thickness  to  chord  ratio  */ 
crp=.6;  /*  fin  percentage  of  total  body  length  */ 
lamda=. 585686;  /*  leading  edge  sweep  angle  */ 

f=l;  /*  nose  fineness  ratio,  (nose  length) / (body  diameter)  */ 

/*  CONSTANTS  */ 

rhoo=1.03;  /*  air  density  at  sea  level  in  kg/m''3  */ 
pi=3. 141593; 

g=9.81;  /*  gravity  in  m/s''2  */ 
gamma=1.4;  /*  gamma  for  air  */ 

R=287;  /*  gas  constant  for  air  in  SI  units*/ 

.  mu=1.8e-5;  /*  coefficient  of  viscosity  */ 
a=340.17; 

/*  Speed  of  sound  at  sea  level  on  standard  day  in  m/s  */ 

/*  INITIAL  VALUES  */ 
x=0.0;  /*  initial  position  */ 
y=0.0;  /*  initial  altitude  */ 
t=0.0;  /*  initial  time  */ 
j=0;  /*  counter  variable  */ 

/*  CALCULATED  INITITIAL  VALUES  AND  CONSTANTS  */ 
vx=v*cos (theta) ;  /*  initial  lateral  velocity  in  m/s  */ 
vy=v*sin (theta) ;  /*  initial  vertical  velocity  in  m/s  */ 

M=v/a;  /*  Mach  Number  for  muzzle  velocity  at  sea  level  */ 

/*  Projectile  Dimensions  and  Parameters  */ 
r=rin*.025417;  /*  radius  in  meters  */ 
r2=r*r; 
r3=r*r*r; 
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stm=stin* . 025417;  /*  skin  thickness  in  meters  */ 
cw=cwin* . 025417;  /*  width  of  contact  area  in  meters  */ 
cw3=cw*cw*cw; 

nl=f*2*rin;  /*  nose  length  in  inches  */ 
nlm=f*2*r;  /*  nose  length  in  m  */ 
nlm2=nlm*nlm; 
nlm3=nlm*nlm*nlm; 

rn=( (nlm*nlm) +r2) / (2*r) ;  /*radius  of  curviture  of  nose  in  m  */ 

rn2=rn*rn; 

rn3=rn*rn*rn; 

thetan=acos (nlm/rn) ;  /*  angle  of  arc  swept  by  nose  */ 
ein=rin-stin;  /*  explosive  radius  in  inches  */ 
er=ein* . 025417;  /*  explosive  radius  in  meters  */ 
f3=f*f*f; 

nv=(pi* ( (nlm* ( (2*rn2) - {2*rn*r) +r2) )-(nlm3/3)-( (rn-r)*nlm* 

sqrt (rn2-nlm2) ) - ( (rn-r) *rn2*asin (nlm/rn) ) ) ) - ( (cw*rn2* 
pi*.5)- (ew*rn2*acos ( (f*cw) / (2*rn) ) ) - ( (2/f) *rn3*sqrt (1- 

( ( {f*cw) / (2*rn) ) * ( (f *cw) / (2*rn) ))))+( (2*rn3) /f ) - ( (f*rn* 
cw/12) * ( (cw/ (2*sin (acos ( (f*cw) / (8*rn) ) ) ) )  +  (cw/ {2*sin 
(acos ( (f*cw) / (4*rn) ) ) ) ) + (3*cw/ {2*sin (acos ( (3*f*cw) / 

(8*rn) ) ) ) )+ (cw/ (2*sin(acos ( (f*cw) / (2*rn) ))))))+(( (cw*f3) / 

(12*rn) ) * ( (cw3/ (128*sin (acos ( (f *cw) / (8*rn) ) ) ) ) + (cw3/ 

(32*sin (acos ( (f*cw) / (4*rn) ))))+( (27*cw3) / (128*sin (acos 
( (3*f*cw) / (8*rn) ) ) ) ) +(cw3/ (8*sin(acos ( (f*cw) / (2*rn) ))))))); 

/*  nose  volume  in  m''3  */ 
nw=rhom*nv;  /*  nose  mass  in  kilograms  */ 

A=(pi*r2)  -  (2*  ( (r2*  (  (pi/2),-acos  ( .  5*cw/r) ) )  -  ( ( (r-  (r*  (1-sin  (acos 
.  ( . 5*cw/r) ) ) ) ) * (r- (r* (1-sin (acos ( . 5*cw/r) )))))*(. 5*cw/r) * 
(1/sin (acos (. 5*cw/r) ))))) ; 

/*  projectile  cross  section  in  m''2  */ 
eGS= (pi* ( (r-stm) * (r-stm) ) ) - (2* ( ( ( (r-stm)  * (r-stm) ) * ( (pi/2) -acos 
( . 5*cw/ (r-stm) )))-((( (r-stm)  - ( (r-stm)  * (1-sin (acos ( . 5*cw/ 
(r-stm) )))))*( (r-stm) - ( (r-stm)  * (1-sin (acos ( .5*cw/ 

(rstm) ))))))* ( .5*cw/ (r-stm) ) * (1/sin (acos ( .5*cw/ (r- 
stm) )))))); 

/*  explosives  cross  section  in  m''2  */ 
el=ew/ (rhoe*ecs) ;  /*  length  of  explosive  in  m  */ 
tlm=nlm+el+stm;  /*  total  length  of  projectile  in  m  */ 
l=(5*r*  (l-sin'(acos  (2/ (5*rin) ) ) )  )  ;  /*  intermediate  calculation  */ 
SA=  (2*'pi*  (  (2*f*r*rn)  -  (rn*  (rn-r)  *asin  (2*f  *r/rn) ) ) )  -  (rn*  ( (pi*cw)  - 
(2*cw*acos ( (f*cw) / (2*rn) ) ) - ( (4*rn*sqrt (1- ( ( (f *cw) / (2*rn) ) 

* ( (f*cw) /(2*rn) )) ) ) /f)+(4*rn/f) ) )+(2* ( (rn2* ( (pi/2) -acos 
(cw/rn)  ) )  -  ( (  (rn-  (cw/2) )  *  (rn-  (cw/2) ) )  *  (cw/rn)  *  (1/ sin  (acos 
(cw/rn) ) ) ) ) )  +  (4*r* (el+stm) tasin ( .5*cw/r) )  +  (2*cw* (el+stm) ) ; 
/*  projectile  surface  area  in  m''2  */ 
sigma=(360*atan(rin/nl) ) /pi; 

/*  semi vertex  angle  at  ogive  point  in  degrees  */ 
br=r- (tlm*. 014054) ; 

/*  base  r  in  m  w/  8  deg  boatail  .1  the  length  of  projectile  */ 
omega= (pi/2) -asin (cw/ (2*r) ) ; 

/*  half  angle  of  base's  curved  section  */ 

ibrh=br*cos (omega) ;  /*  intermediate  step  */ 
ibrw=br*cos ( (pi/2) -omega) ;  /*  intermediate  step  */ 

Sb= (2* ( (br*br* ( (pi/2) -omega) ) - (ibrh*ibrh* (cos (omega) 

/sin(omega) ) )  ) ) +  (2*cw*ibrw)  ;  /*  base  area  in  m''2  */ 
cr=crp*tlm;  /*  root  chord  in  meters  */ 

psi= (pi/2) -lamda;  /*  "inside"  sweep  angle  of  leading  edge  */ 
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b=2* ( (sqrt (2* ( {r*sin(acos ( .5*cw/r) ) ) * (r*sin(acos 

{.5*cw/r) ) ) ) ) )-r) ;  /*  span  of  single  fin  in  meters  */ 

ct=(cr-(.5*b/tan(psi) ) ) ;  /*  tip  chord  in  meters  */ 

Psi=pi-psi;  /*  "inside  sweep’  angle  of  tip  */ 

sf=  { .  5*b*  (cr-ct) )  +  (ct*b) ;  /*  planform  area  of  fin  pair  in  m''2  */ 

lamdam=ct/cr;  /*  intermediate  calculation  *! 

mac=  ( .  666667*cr*  (l+lamdam+  {lamdam*lamdam) ) )  /  (1+lamdam) ; 

/*  mean  aerodynamic  chord  in  meters  */ 
tmac=miac*tau;  /*  fin  thickness  at  MAC  in  meters  */ 
rf=mac*(l+(tau*tau))/(4*tau) ; 

/*  radius  of  curviture  of  biconvex  fins  at  MAC  in  meters  */ 
phi= {pi/2) -asin (mac/ (2*rf) ) ; 

/*  lower  limit  of  integration  for  quarter  of  fin  volume  */ 
fv=2*b* { ( (rf*rf ) * ( (pi/2) -phi) )- ( ( (rf*rf) - (rf*tmac) + ( .25*tmac 
*tmac) ) *cos (phi) /sin (phi) ) ) ; 

/*  volume  of  single  fin  in  m''3  */ 
fw=2*rhom*fv;  /*  mass  of  all  four  fins  in  kg  */ 
fsa=16*rf* ( (pi/2) -phi) *b;  /*  surface  area  of  four  fins  */ 
m= (rhom* ( ( (A-ecs) *el)  +  (stm*A) ) ) +ew+nw+fw 
/*  projectile  mass  in  kg  */ 
dt=.0001* (v/g) ;  /*  time  step  for  integration  */ 

printf ("%f , %f ,  %f , %f , If \n", x/1000, y/1000, v, t,M) ; 
while{y>=0  &&  y<11000) { 

T=288- ( (6.5*y) /lOOO) ;/*  Adiabatic  Cooling  w/  Alt  Increase*/ 

a=sqrt (gamma*R*T) ; 

dc= { {y*g) / (R* (T-288 . 16) ) ) +1; 

dce=pow( {T/288.16) ,-dc) ; 

rho=rhoo*dce; 

v=sqrt (vx*vx+vy*vy) ; 

M=v/a; 

beta=sqrt( (M*M)-1) ; 

Mu=asin(l/M) ; 

Dl=tan (psi-Mu) ;D2=tan (Psi-Mu) ;D3=D1+D2; 

/*  intermediate  calcs  */ 

Reb=(rho*v*tlm) /mu;  /*  Reynolds  nvimber  for  body  */ 

Ref=(rho  *v*mac) /mu;  /*  Reynolds  number  for  fins  */ 
pb=pow( (log(Reb) /log{10) ) ,-2.58)  ; 
pf=pow( (log(Ref ) /log(lO) ) ,-2.58) ; 

Cfb=-455*pb;  /*Yellow  book*/ 

Cff=.455*pf; 

comf=pow( (1+ ( ( (gamma-1) /2) *M*M) ) ,  - . 467) ; 

/*  coitpressible  flow  factor  */ 

CDfb= (Cfb*SA*comf ) /A; 

CDf f= (Cf f *fsa*comf ) /A; 
sigmac=pow{ (sigma/10) ,  1.69) ; 
pr=sigmac* ( . 083+ { . 096/ {M*M) ) ) ; 

CDFb=pr* (1- ( ( (196*f*f ) -16) / (14* (M+18) * (f*f ) ) ) )  ; 
seff=sf-( ( (1.1452*D1/D3)*(1.1452*D1/D3)*D2)+ 

( (1. 1452*D2/D3) * (1.1452*D2/D3) *D1) ) ; 
CDFf=(5.33*tau*tau* (sf- ( .5*seff ) ) ) / (beta*sf ) ; 
lbt=500000* ( ( . 0555*M*M)  +  ( . 006*M)  +  (- . 00255*M*M*M) +. 85) ; 
nbt= ( (-485*M) * (-485*M)  ) - (Ibt) ; 

CDbt=(10* { {485*M)-sqrt (nbt) ) /250000) -.00375; 
lb=4000* ( { . 007*M*M)  +  { . 00007*M)  +  (- . 0009*M*M*M) +. 85) ; 
nb={ {-445*M)*(-445*M) )-{lb) ;  ’ 

CDb=110*Sb* ( (445*M)  -sqrt (nb) ) / (2000*A) ; 
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C=CDf b+CDf f +CDFb+CDFf +CDbt+CDb ; 
k=(C*rho*A) / (2*m) ; 
vx=vx- (k*v*vx*dt) ; 
vy=vy-g*dt- (k*v*vy*dt)  ; 
x=x+vx*dt ; y=y+vy*dt ; t=t+dt ; 
j=j+l; 
if(j>=100) { 

printf ("%f,%f,%f,%f,%f\n",x/1000,y/1000,v,t,M) ; 
j=0;}} 

while (y>=11000  &&  y<25000) { 

Isothermal  Region 
Equations  of  Motion  Loop} 
while (y>=25000  &&  y<47000) { 

T=216. 5+ (3* (y-25000) /lOOO) ;  /*  Adiab.  Heating  w/  Alt  Inc* 
Equations  of  Motion  Loop} 
while (y>=47000  &&  y<53000) { 

Isothermal  Region 
Equations  of  Motion  Loop} 
while {y>=53000  &&  y<79000) { 

1=282.5- (4 .5* (y-53000) /lOOO) ;  /*  Adiab.  Cooling  w/  Alt 

Inc*/ 

Equations  of  Motion  Loop} 
while (y>=7 9000  &&  y<=90000) { 

Isothermal  Region 
Equations  of  Motion  Loop} 
while (y>=90000  &&  y<=105000) { 

T=165.5+(4* (y-90000) /lOQO) ;  /*  Adiab.  Heating  w/  Alt  Inc*/ 
Equations  of  Motion  Loop} 

while (y>=105b00) {Drag  is  Effectively  Zero  Above  this  Altitude  so 
all  the  drag  calculations  are  removed  from  the  loop  and  the 
drag  terms  are  removed  from  the  equations  of  motion} 
while (y>=90000) { 

Isothermal  Region 

T=165.5+{4* (y-90000) /lOOO) ;  /*  Adiab.  Cooling  w/  Alt  Dec*/ 
Equations  of  Motion  Loop} 
while (y>=79000) { 

Isothermal  Region 
Equations  of  Motion  Loop} 

'  while (y>=53000) { 

T=282.5- (4.5* (y-53000) /lOOO) ;/*  Adiab.  Heating  w/  Alt  Dec*/ 
Equations  of  Motion  Loop} 
while (y>=47000) { 

Isothermal  Region 
Equations  of  Motion  Loop} 
while (y>=25000) { 

T=216.5+ (3* {y-25000) /lOOO) ;  /*  Adiab.  Cooling  w/  Alt  Dec*/ 
Equations  of  Motion  Loop} 
while (y>=11000) { 

Isothermal  Region 
Equations  of  Motion  Loop} 
while (y>=0) { 

T=288- { (6.5*y) /lOOO) ;  /*  Adiabatic  Heating  w/  Alt  Dec*/ 
Equations  of  Motion  Loop} 
xf=x- ( {vx*y) /vy) ; 
yf=y- ( (vy*y) /vy) ; 

printf  (”%f,%f,%f,%f,%f\n’*,xf/1000,yf/1000,v,t,M) ;} 
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APPENDIX  (C)  MATLAB  STABILITY  PROGRAM 

%  stability  Study 
%  Inputs 

r=2-5;  %  projectile  radius  in  inches 
h=2;  %  width  of  contact  area  in  inches 
f=l;  %  nose  fineness  ratio 
st=.5;  %  skin  thickness  in  inches 

rhos“19200.0;  %  density  of  metal  in  kg/m''3  (tungsten) 

rhoe=1520.0;  %  density  of  explosive  in  kg/m‘'3  (C4) 

ew=5.5;  %  explosive  weight  in  kg 

tau=0.03;  %  fin  thickness  ratio 
z=.6;  %  root  chord/total  length* 

%  Constants 
gainma=l‘.  40; 

rho=. 00004;  %  air  density  in  Ibm/cubic  inch 

T=520;  %  Standard  ten¥)erture  at  sea  level  in  Rankine 

a=(sqrt (gamina*32.2*53.3*T) ) *12;  %  Speed  of  sound  in  inches/second 

%  Nose/Body 

d=(2*r);  %  Projectile  width  in  inches 

rm=r*. 0254166667;  %  radius  in  meters 

hm=h*. 0254166667;  %  width  of  contact  areas  in  meters 

stm=st*. 0254166667;  %  skin  thickness  in  meters 

nl=2*rm*f;  %  nose  length  in  meters 

nli=2*r*f;  %  nose  length  in  inches 

Rm=(  (nl''2)  +  (rm''2) )  /(2*rm)  ; 

%  Radius  of  cuirvature  of  nose  in  meters 
thetas=acos (nl/Rm) ; 

nv=  (pi*  ( (nl*  ( (2*  (Rm^2) )  -  (2*Rm*rm)  +  (rm''2) ) )  -  ( (nl^3)  /3)  -  (nl*  (Rm-rm)  * . . 
sqrt  ( (Rm''2)  -  (nl''2) ) )  -  { (Rm-rm)  *  (Rm.''2)  *asin  (nl/Rm) ) ) )  -  ( (hm*  (Rm'^2)  * 
pi*.  5)  -  (hm*  (Rm^2)  *acos  ( (f*hm)  /  (2*Rm) ) )  -  ( (2/f )  *  (Rm^3)  *  (sqrt  (1-. . . 

(  ( (f*hm)  /  (2*Rm) )  ^^2) ) ) )  +  ( (2*  (Rm''3) )  /f )  -  ( (f*Rm*hm/12)  *  ( (hm/  (2*sin. . 
(acos ( (f*hm) / (8*Rm) ) ) ) )  +  (hm/  (2*sin (acos ( (f *hm) / (4*Rm) ) ) ) )  +  (3*hm. / 

(2*sin(acos( (3*f*hm) / (8*Rm) ) ) ) )  +  (hm/(2*sin(acos( (f*hm) /(2*Rm) ))))))  + 
( ( (hm*  (£"^3) )  /  (12*Rm) )  *  ( ( (hm''3)  /  (128*sin  (acos  ( (f*hm)  /  (8*Rm) ))))  +  .. 
( (hm''3)  /  (32*sin(acos  ( (f  *hm)  /  (4*Rm) ))))  +  ( (27*  (hm'^S) )  /  (128*sin  (acos 
( (3*f*hm)  / (8*Rm) ))))  +  ( (hm^3) /  (8*sin(acos ( (f*hm) / (2*Rm) ))))))); 

%  Nose  volume  in  m'^S 

nw=rhos.*nv;  %  Weight  of  nose  in  kilograms 

pcs=(pi*  (rm''2)  )-(2*(  ((rm''2)*(  (pi/2) -acos  (. 5*hm/rm) ) )-( ( (rm-(rm*(l- 

sin. . . 

(acos  ( .5*hm/rm) ) ) ) )  ^^2)  *cot  (acos  ( .5*hm/rm) ) ) ) )  ; 

%  Projectile  cross  section  in  meters''2 
pcsi=pcs/. 000646;  %  Projectile  cross  section  in  in^2 
ecs=(pi*(  (rm-stm)  . '>2) )  -  (2*  ( ( ( (rm-stm) '‘2)  *  ( (pi/2) -acos  ( .5*hm/ (rm- 
stm) )))-... 

( (  ( (rm-stm)  -  ( (rm-stm)  *  (1-sin  (acos  ( .  5*hm/  (rm-stm) ) ) ) ) )  "^2)  *cot  (acos 
(. 5*hm/ (rm-stm) ))))) ;  %  Explosives  cross  section  in  m^2 

%  Explosives  cross  section  in  m''2 
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el=ew/ (rhoe*ecs ) ;  %  Length  of  explosive  in  meters 
eli=el/ . 025416667;  %  Length  of  explosive  in  inches 
pwk=(rhos*( ( (pcs-ecs) *el)+(stm*pcs) ) )+ew+nw; 

%  Projectile  weight  in  kilograms 
pwl=pwk*2.2;  %  Projectile  weight  in  Ibm 

tlm= (f *2*rm) +el+stm;  %  Total  length  of  projectile  in  meters 
tli=tlm/ . 0254166667;  %  Total  length  in  inches 

sam=  (2*pi*  (  (2*f  *rm*Rm)  -  (Rm*  (Rm-rm)  *asin  (2*f  *rm/Rm) ) ) )  -  (Rm*  ( {pi*hm)  - 

(2*... 

hm*acos  (  (f  *hm)  /  (2*Rm)  )  )  -  (  (4*Rm*sqrt  (1-  (  ( {f*hm)  /  {2*Rm) )  ^2) )  )  /f )  +. . . 
(4*Rm/f ) ) )  +  (2*  ( ( (Rm"2)  *  (  (pi/2)  -acos  (hm/Rm) ) )  -  { (  (Rm-  (hm/2)  )  '^2}  *cot . . . 
(acos  (hm/Rm)  ) )  ) )  +  (4*rm*  (el+stm)  *asin  ( .5*hm/rm)  )  +  (2*hm*  (el+stm) )  ; 

%  Projectile  Surface  Area  in  m^2 
sai=sam. /. 000646;  %  Projectile  Surface  Area  in  in'^2 

%fins 

b=2*  ( (sqrt  (2*  ( (r*sin  (acos  ( .  5*h/r) ) )  ^^2) ) ) -r) ;  %  span  in  inches  w/o  body 
bm=b* . 025416667 ;  %  span  in  meters 
cr=2*tli;  %  centerline  chord  in  inches 
lamda=. 585686; 

%  sweep  ang.  of  34  deg (in  rads);  keeps  LE  supersonic  down  to  Mach  1 
theta=  (pi/2) -lamda;  %  *’inside"  sweep  angle  of  56  degrees 
ct= (cr- (. 5*b/tan (theta) )) ;  %  tip  chord  in  inches 

Sf= ( . 5*b* (cr-ct) ) + (ct*b) ;  %  planform  area  of  fin  pair  in  square 

inches 

Sd= (cr^2) *tan (theta) ;  %  area  if  delta  where  whole  in  sgare  inches 

Sl=Sf/Sd; 

S2=Sf /pcsi; 

mac=(2*cr*(H-(ct/cr)  +  (  (ct/cr)  "^2)  )  )  /  (3*  (1+ (ct/cr)  )  )  ; 

%  mean  aerodynamic  chord  of  fins  in  inches 

macm=mac* . 025416667;  %  mean  aerodynamic  chord  of  fins  in  meters 
Rmf=(  ( ( .5*mac)  ^2)  +  (  ( .  5*tau*inac)  ^2) )  /  (tau*mac)  ; 

%  radius  of  curviture  of  biconvex  fin  at  mac  in  inches 
Rmfm=Rmf*. 0254166667 ; 

%  radius  of  curviture  of  biconvex  fin  at  mac  in  meters 
thetaf=acos ( . 5*mac/Rmf ) ; 

fw=2^rhos*bm*  ( (  (Rmfm"^2)  *  (  (pi/2)  -thetaf )  )  -  ( ( (Rmfm-  ( .  5*tau*macm)  )  "2)  .  .  . 

*cot (thetaf ) ) ) ;  %  Weight  of  single  fin  in  kg 
phi=pi-theta;  %  "inside”  tip  angle  of  129  degrees 
AR-b'^2/Sf; 

Kl-5.33; 

Apri=2*tau/3; 

%  Center  of  Gravity 

tw=pwk+ (2*fw) ;  %  weight  of  entire  projectile  including  fins  in  kg 
CGn=.707*nli*nw;  %  CG  point  for  the  nose 

CGb={ ( .5*eli)+nli) * ( (rhos* (pcs-ecs) *el)+ew) ;  %  CG  point  for  body 
CGep=( (•5*st)+eli+nli) * (stm*pcs*rhos) ;  %  CG  point  for  the  end  plate 
CGf-( (mac*.5)+( (l-(z+,l) )*tli))*2*fw; 

CG= (CGn+CGb+CGep+CGf ) /tw;  %  CG  point  of  projectile 

%  Interference  Factors 
rs=r/ (r+(b/2) ) ; 

Kbf=(.869*(  (rs+.65)"^2)  )-.367; 

Kfb=(.278*(  (rs+1.3)"^2)  )+.529; 
kbf=rs; 

kfb-{.2*((rs-.5)-2))+.95; 
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M=linspace (1.3, 9, 250) ; 
beta=sqrt  (M.  ''2-1) ; 
mu=asin ( 1 . /M) ; 
if  (beta./f)<=l 

CNan=(-.28.* { ( (beta. /f) -1.9) .^2) )+3.58; 

%  Approximation  over  Mach  range  for  chart  pg  13-4 
Xn=CG- (tli.* ( (.125.*beta./f)+.15) ) ; 

%  From  graph  pg  14-8  adjusted  to  CG 
else 

CNan=(-3.4.*( ( (f . /beta) - . 9) .^2))+3.4; 

Xn=CG-(tli.*(  (-.125'.  *f. /beta) +.425) ) ; 
end 

Dl=tan (theta-mu) ; 

D2=tan (phi-mu) ; 

D3=D1+D2;- 

if  ( (b/2) - ( (cr+ct) .*tan(mu) ) )<=0 

Seff=Sf-( ( ( (1.1452.*D1./D3) ."2) .*D2)+( ( (1.1452.*D2./D3) .^2) .*D1) ) ; 
%  Effected  tip  area  (split  into  two  terms  to  account  for  overlap) 

Seff=(  (cr''2)  .*tan(mu) )  +  (  (ct^2)  .*tan(mu) )  ; 

end 

CNa2df=(4./beta) ; 

%  CNalpha  before  tip  effects 
CNaf=CNa2df-( (2.*Seff) ./(beta.*Sf) ) ; 

CNa=CNan+ (CNaf . * (Kbf+Kfb) . *S2) ; . 

CNd=CNaf  .*  (.(kbf+kfb)  *S2)  ; 

Xf=-(cr*. 666666* (SI- ( . 5. *Seff . /Sd) ) ) - ( ( (1- (z+. 1) ) *tli)-CG) ; 

%  CP  of  fin  adjusted  for  missing  area  and  to  CG 
Xbf=-(cr.* ( (.9g88;*beta.*AR)+.35) )-( ( (l-(z+.l) ) *tli)-CG) ; 

%  From  NACA  figure  15i  and  adjusted  to  CG 

CMa= ( (CNan. *Xn) + (CNaf . * ( ( (Xbf . *Kbf ) - (Xf . *Kfb) ) . *S2) ) ) . /d; 

CMl=CMa .*.017453; 

CM3=CMa.*. 052360; 

CM5=CMa.*.087227; 

CM10=CMa.*.174533; 

CMd=(CNaf.*((Xbf.*kbf)+(Xf*kfb) )*S2) ./d; '%From  pg  20-2 
q=.5.*rho.*((a.*M).'^2); 

CN=pwl./(q.*pcsi) ; 

atrim= ( (CN . *CMd) . / ( (CNa . *CMd) + (CNd. *CMa) ) ) . *57 . 3; 

LFS=(CNd)  .*q.*pcsi./pwl; 

cpn=Xn . /d; 

cpf=Xf./d; 


figure ( 1 ) , plot (M, CNa, ' k ' ) , grid 

title ('CN  alpha '),xlabel(' Mach '),ylabel(’CN  alpha  (/rad)') 
figure ( 2 ) , plot (M, CNd, ' k ' ) , grid 

title ( ' CN  delta ') , xlabel ( ' Mach '), ylabel (' CN  delta  ( /rad) ' ) 
figure(3),plot(M,CMa, 'k'),grid 

title (' CM  alpha ' ) , xlabel ( ' Mach ' ) , ylabel ( ' CM  alpha  ( /rad) ' ) 
figure (4) ,plot (M, CMd, ' k' ) , grid 

title ('CM  delta'), xlabel ('Mach'), ylabel (’CM  delta  (/rad)') 
figure ( 5 ) , plot (M, atrim, ' k ' ) , gr id 

title ( 'Alpha  trim' ), xlabel ( 'Mach' ) , ylabel ( 'Alpha  trim  (deg) ' ) 
figure  (6)., plot  (M, LFS, '  k' )  ,grid 
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title ( *  Load  Factor  Sensitivity ’ ) , xlabel ( *  Mach ’ ) , ylabel ( ’ (n~l ) /delta ’ ) 
figure(7) ,plot (M,CNan, ’k’) ,grid 

title < ’Cnalpha  Nose/Body  Alone ’) f xlabel (* Mach M / ylabel (’ CNalpha 
(/rad) ’) 

figure (8) ,plot (M, CNaf ^  ) , grid 

title { ’Cnalpha  Fin  Alone ’ ) ^xlabel ( ’Mach’ ) , ylabel ( ’ CNalpha  (/rad) ’ ) 
figure  ( 9 )  plot  (M,  cpn,  ’  k  ’ ) ,  grid 

title (’Xcp  Nose/Body ’), xlabel ( ’Mach’ ), ylabel { ’Xcp (from  eg) /body 
diameter’ ) 

figure (10) ,plot (M, cpf , ’ k’ ) , grid 

title ('Xcp  Fin’) ^xlabel ('Mach’ ), ylabel (’Xcp (from  eg) /body  diameter’) 
figure (11) ^plot (M,CMl, 'k' ) /grid 

title (’CM  at  One  Degree’ ) /Xlabel ( 'Mach' ), ylabel ( 'CM' ) 
figure ( 12 ) / plot (M/ CM3 / ' k ' ) , grid 

title (’CM  at  Three  Degrees ’)/ xlabel ( ’Mach' ), ylabel ( ’CM’ ) 
figure (13) /plot (M, CM5^  *  k ’ ) , grid 

title (’CM  at  Five  Degrees ’)/ xlabel ( ’Mach’ )/ ylabel ( ’CM’ ) 
figure (14) /plot (M, CMIO/ ’ k’ ) / grid 

title (’CM  at  Ten  Degrees ')/ xlabel ( ’Mach’ )/ ylabel (’ CM’ ) 
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APPENDIX  (D)  Graph  Approximation  Pr<^rams 


The  first  program  provided  approximates  the  four  interference  fectors. 

%  Interference  Factors 


clear 

rs=linspace (0,1,20); 


%KBW 

A=.869; 

B=.65; 

C=-.367; 

KBW=(A.* { (rs+B) .^2) )+C; 

%kBW 
kBW=rs ; 

%KWB 

A=.278; 

B=1.3; 

C=.529; 

KWB={A.*( (rs+B) .^2) )+C; 


%kWB 

A=.2; 

B=-.5; 

C=.95; 


kWB=(A.*  ( (rs+B)  .''2)  )+C; 
figure (1) 

plot(rs,KBW,  'k', .2,  .279,  'ok',  .4, .607,  ’ok’, .6,1, ’ok', 
rs,kBW, 'k’, .2, .213, ' dk’ ,  . 4, . 41,  ' dk’ , . 6, . 61, ' dk’ , . 
rs,KWB, ’k’,  .2,1.16,  ’*k’,  .4,1.35,  ’*k’, .6,1.55, ’*k' 


.8,1.446, ’ok’,  . . . 

8, .8, ’dk’,  . . . 

,  .8,1.77, ’*k’),grid 


hold 

plot(rs,kWB,  ’k’,  .2,  .958,  ’^'k’,  .4,  .933,  ’''k' 

title ( ’ Interference  Factors ’ ) 

xlabel (’ Radius  to  Semispan  Ratio’) 

text(.43,1.5, ’KWB’) 

text (.64, 1.3, ’KBW’) 

text(.21, .91, ’kWB’) 

text(.41, .35, ’kBW’) 


6,  .94,  ’^k',  .8,  .967,  ’'"k’) 


Figiire(19)  on  the  next  page  is  the  plot  of  this  program 
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Figure  (19)  interierence  Factors 
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The  nose/body  lift  coefficient  is  calculated  by  two  programs.  The  curve  changes 
when  the  the  fineness  ratio  equals  p  and  the  two  branches  needed  to  be  done  seperately. 

%  Nose-Body  CNalpha  Graphic  Approximation 

%  Left  side  of  USAF  DATCOM  figure  4.2.1.1-5a  on  pg  13-4  of  AA3701  notes 

A=-.28; 

B=-1.9; 

C=3 . 58 ; 

Beta=linspace (.4,1, 100) ; 

CNalpha=(A.*(  (Beta+B)  .''2)  )+C; 
figure (1) 

plot (Beta, CNalpha, ' k’ , . 4, 2 . 95, ' ko’ , . 6, 3. 1, ’ ko' , . 8, 3 .25, * ko',!, 3 . 35, 'ko' 
) , grid 

title {' CNalpha  vs  Beta/fN') 

xlabeK'Beta  divided  by  Nose  Length/Body  Diameter  ratio') 
ylabel ( ' Cnalpha ' ) 


%  Nose-Body  CNalpha  Graphic  Approximation 

%  Right  side  of  DSAF  DATCOM  figure  4.2.1.1-5a  on  pg  13-4  of  AA3701 
notes 

A=-3.4; 

B=-.9; 

C=3.4; 

Beta=linspace ( . 1118, 1, 100) ; 

CNalpha= (A. * ( (Beta+B) . ^2) ) +C; 
figure (1) 

plot (Beta, CNalpha, 'k', .5,2.86, 'ko', .6,3.09, 'ko', .8,3.36,  'ko',  .9,3.4,  ko 

t 

1,3.35, 'ko' ) ,grid 
title ( 'CNalpha  vs  Beta/fN') 

xlabel('Beta  divided  by  Nose  Length/Body  Diameter  ratio') 
ylabel ( ' Cnalpha ' ) 


Figure  (20)  and  Figure  (21)  on  the  next  two  pages  are  the  plots  of  these  two  programs. 
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Cnafpha 
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Cnalpha 


79 


The  boattail  drag  and  base  drag  coefficients  are  calculated  from  the  following 


graphs. 

A=-2000; 

B=65000; 

C=.69; 

D=-570; 

E=l; 

F=.0115; 

G=-.00097; 

M=linspace ( 1 . 2 , 9 , 100 ) ; 

ra=(  {D.*M)  . ''2)  -  (4  .  *B.  *  ( (C.  *M)  +  (F.  *  (M.  ^3) )  +  (G.  *  (M.  ^4)  ) +E) )  ; 
cd={10.* (-(D.*M)-sqrt (m)  ) ./ (2*B) )-.015; 
figure (1) 

plot  (M,cd,  'kM.5,  .0171,  '  ko' , 2,  .  01,  '  ko' ,  2 . 5,  .0071,  ’  ko' ,  3, . 0057,  '  ko’ ,  . . 

3.5, .0054, 'ko',4, .0051, ’ko*), grid 
title ( 'Boattail  Drag  Coefficient  vs  Mach  Number') 
xlabelCMach  Number '), ylabel  (' Boattail  Drag  Coefficient') 


%  Base  Pressure  Coefficient  Graphic  Approximation 

%  From  Chin  pg  31 

A=.007; 

B=1000; 

C=. 00007, • 

D=-445; 

E=.85; 

F=-.0009;  •  • 

M=linspace (1, 9, 100) ; 

m=(  (D.*M)  .''2)-{4.*B.*  (  {A.*(M.''2)  )  +  (C .  *M)  +  (F.  *  (M.  ^'S)  ) +E)  )  ; 
cd=110.* (- (D.*M) -sqrt(m) ) ./{2*B) ; 
figure (1) 

plot(M,cd, 'k',1.2, .168, 'ko',1.4, .148, 'ko',1.8, .12, 'ko',2.2, .1, 'ko',2.6 
.082,  . . . 

.  'ko', 3.0, .072, 'ko'), grid 

title ('Base  Pressure  Coefficient  vs  Mach  Number') 
xlabel ( 'Mach' Number' ), ylabel ( 'Base  Pressure  Coefficient') 


The  plots  are  presented  as  Figure  (22)  and  Figure  (23)  on  the  next  two  pages. 
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BtmtUill  Drao  Coefficient 


81 


LIST  OF  REFERENCES 


1 .  Chin,  S.S.,  Missile  Configuration  Design,  McGraw-Hill  Book  Conq)any,  Inc.,  1961. 

2.  Pitts,  W.  C.,  Nielsen,  J.  N.,  and  Kaattari,  G.  E.,  Lift  and  Center  of  Pressure  of  Wing- 
Body-Tail  Combinations  at  Subsonic,  Transsonic,  and  Supersonic  Speeds,  NACA 
Report  1307. 

3.  Kemp,  N.  H.  and  Riddell,  F.  R.,  "Heat  Transfer  to  Satellite  Vehicles  Re-entering  the 
Atmosphere,"  Jet  Propulsion,  vol.  27,  no.  2,  part  1,  p.  137,  February  1957. 

4.  Zucker,  R  D.,  Fundamentals  of  Gas  Dynamics,  Matrix  Publishers,  Inc.,  1977. 

5.  Hoak,  D.  E.  and  Finck,  R.  D.,  U.S.  Air  Force  Stability  and  Control  DATCOM  Vol3 
5ec.  5,  Engineering  Documents,  1960. 

6.  Blackadar,  A.  BC.,  "Wind,"  Microsoft  Encarta  97  Encyclopedia,  Microsoft 
Corporation,  1996. 

7.  Anderson,  J.  D.  Jr.,  Introduction  to  Flight,  McGraw-Hill  Book  Company,  Inc.,  1979. 


83 


84 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center  .  2 

8725  John  J.  Kingman  Road,  STE  0944 

Ft  Belvoir,  V  irginia  22060-62 1 8 

2.  Dudley  Knox  Litoary . . * .  2 

Naval  Postgradu2de  School 

411  Dyer  Rd. 

Monterey,  California  93943-5 101 

3.  Engineering  &  Technology  Curricular  Office  (Code  34)  .  1 

Naval  Postgraduate  School 

Code  34 

700  Dyer  Road,  Bldg  245  Rm  1 15 
Monterey,  California  93943-5109 

4.  Professor  Conrad  F.  Newberry  .  2 

AA/Ne 

Department  of  Aeronautics  and  Astronautics 
Naval  Postgraduate  School 
Monterey,  California  93943-5106 

5.  Professor  William  B  Maier,  n  .  2 

Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  California  93943 

6.  LCDR  Francisco  M.  Juanche .  1 

1050  Camino  Atajo 

Chula  Vista,  California 


85 


